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INDEFINITE QUADRATIC FORMS IN MANY VARIABLES 


Mary E. FLAHIvE 
Department of Mathematics, University of Lowell, Lowell, Massachusetts 01854 


(Received 15 December 1983; after final revision 6 June 1988) 


In this paper we investigate the positive and the nonnegative inhomogeneous 
spectra for indefinite binary quadratic forms of rank at least 21. (Recent 
work of Margulis should allow the restriction m A 21 to be reduced.) For 
the positive spectra we prove that the sharpest constant depends only on the 
congruence class (mod 8) of the signature of the form; all extreme forms and 
grids are also obtained. In addition we bound the best possible constant for 
all nonnegative inhomogeneous spectra. 


1. INTRODUCTION 
Let f(X) = x= = ay X1 X, be a nonsingular quadratic form with real coeffi- 
1<i<j<n 
cients. 


Using the notation in Watson”, we let 


2441 Ais, -s>s, “Gin 
Ate VP iigt Yl ites) Jz Oo, 
Gin Aan ee6 24nn 


be called the matrix associated with f and set D = det f to equal the determinant of A. 
Also, the discriminant of f is defined to be 


det f if n is even 
d(f) =(—1)""1. . det f 


> mat | 
D if n is odd (1) 


In Watson?! (p. 3), it is shown that the discriminant of any integral quadratic 
form is an integer; the nonsingularity of fimplies that each of det fandd(/) is 
nonzero. Further, in Watson?! (p. 3), it is shown that, if f, +f, is a disjoint sum, 
then 


d(fi)d(f2) if m nz is even 

Bat Jy) = eye (fi) d (f2) ifm, nz is odd. 
Also, d(/) is an integer and by Watson*? (p. 21), 

d(f) =0, 1 (mod 4), if n is even. t3) 


ky 
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For any V € R", we define 
P+ (f, V) = inf {f(X + V) > 0: X is an integral point}; 
P(f, V) = inf {f(X + V) > 0: X is an integral point} 
where for the special case V = 0 we restrict to X 4 0. We define 
Pr “(f) == sup ten; Vy nv-e RY: 
Fu (f) = sup{P (fV)2 V ER} 


which may be called, respectively, the positive and nonnegative nonhomogeneous 
minima of the form/. For fixed n > 1, | s | <n, we consider the sets 


{P; *(f)'] | D| : findefinite of rank n and signature s}; 
{P: (f)"/ | D | : f indefinite of rank n and signature s} 
which we shall call, respectively, the positive and nonnegative inhomogeneous spectra. 


Bambah et al.** have considered the positive spectrum. They show that, for 
n> 2ands=0,1,2,3:n >7ands— — fs 


Pe ( S| DVS OL sila yet 


They also prove that this bound is best possible for each of these values of n, s. Hence, 


this extends the work of many authors, among them : Davenport and Heilbronn”, 
Blaney®, Barnes‘, Dumir!4"!>, Hans-Gill and Raka!8. 


In addition, Dumir and Hans-Gill'® showed that | is the best constant forn = 4, 


Ss = — 2; in Bambah et al.*, (7/8)° is proved to be best possible forn = 5,5 = — }, 
In this paper we investigate these spectra form > 21. We prove that, for any 

n, S, the best possible constants depend only on s (mod 8); for the various nonnegative 

spectra an upper bound of 1/5 is obtained. The statements of these resul 


ts are given 
in Theorems | to 3 below. 


For convenience, we set 


E, = x}; Ey = x3 + x, x. + x?: 


Es = xi + x? + x3 + XX, + X, Xs + X, Xs; 


Ey = xi + x2 + x? + x? + x§ LET hig 7 ig ib.) Ye + (x, mrt eee 


Ey oat E; (x,, X25 X3, X4, 0); 


iy 
l 


4 4 
+ ; a (xi2 + (2Xi44 + xX — : x x;)); 
—_ =1 


E, = E, (Xie. X55 X65 X6, Xe); 


a 
l 


Ey (0, 05% sy hare 
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These forms are our building blocks; each E; is a positive definite integral quad- 














ratic form with minimum nonzero value equal to !. Also, 
TABLE I 
7 
f E, E; E; E; Es fos Es Es 
det f 2 3 4 4 4 3 4 1 
d(f) 1 —3 —2 4 2 -3 -4 1 





Define each of F, and H, to be the form which is identically zero, and for any 
nonzero integer m, define 


Fin = >> Es (Xakcrsetsh tatty 
0<ixm 


for any positive integer ¢, let 


H, = >> Xoi-4 Xai. 
1<SiSt 


Then 
s (F,,) = 8m, det Fm = 1; s (Hi) = 0, det He = (— 1). 


Definition—f is said to be (integrally) equivalent to g [denoted by f ~ g] if there 
exists an integral matrix N, with det N = + 1, such that f(N X) = g (X). 


Theorem aes be an indefinite quadratic form of rank n > 21 and signature 
s=8&tk,-3<k <4. Then for all V € R*® 


P*(f, V)"]| D| < 1Jmin{ | k| + 1, 4}. 


Moreover, the equality sign is required only when either 


(i) fis: equivalent oe a positive multiple of sgn (7) - Fm + H: + g, for 21 =n 


— &m—n (2) and m and ey are given in Table II; under this equivalence V becomes 
O (mod !) 


TABLE II 
k 0 axl! +2 gs +3 +3 4 
s all all be~ = n—2 |s|=an-2 |s' <n-4 |s|> n—-4 all 
m 1q| 1q| lq| lo) -1 Iq | eas Se ee a 
g 0 +E, +Ee FE6 +E; FE; E, 


Be 


or 
(ii) f is equivalent to a positive multiple of — Heth Pala g x ee. +X, +, +h, 
for suitable choice of ft; and ¢,, and V becomes (1/2, V:) (mod |), where h and Ve are 


given below (Table III) : 
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TABLE III 
k 0 +1 +2 a 4 
h 0 +2F, +F2E; +2E, -2E¢ +4E, +2E,; F2E; 4H, 2E, 2£E¢ 
V2 (mod 1) ... 0 \ 0 0 + 0 0 (7) 0 $(0,1,1,1) 





This extends the following theorem proved by Watson?®. [Note: The fourth excep- 
tional form was not recognized until Watson**]. 


Let f be an integral, primitive, nonsingular quadratic form of rank n and signa- 
tures (f) = 8m+k,for—3<k <4. If f is definite, we assume it is positive 
definite and of rank n < 8. Then, if fis equivalent to none of the forms 4x, x, — x3, 
16x, Xo —xj, aS Xo a xs rae 3X3 X4 — Gx7; Xo (6x; a; Xo te 3x3) — a 

P* (f, 0) n ie: 
Dut ~~ min{|k| + 1, 4° 
The sign of equality is necessary only if either 
| det f| = min {| k| + 1, 4} orf ~ 8x, x, — x2, 


In Watson**, the author finds other elements of the spectrum {Pt (f, 0)"/| D|: 
n(f) =n}, forn = 3. In thatarticle, he also gives an inductive argument for esti- 


mating P* (f,0)"/ | D| for zero forms of rank n > 4 from information on forms of 
rank n — 2. 


Theore m 2—Let f be an indefinite quadratic form of rank n > 21 and signature 
s=8¢ + k, —-3<k <4. If fis none of the exceptional forms listed in Theorem I, 
then P* (f, V)"] | D| < 1/4. Moreover, for|k] <2 equality is required only if 
either 

(i) f is equivalent to a positive multiple of sgn (q)- Fm + Hy; + 5 for 2t = n 
— 8m — n(g) and m and g are given in Table IV. 

TABLE IV 
0 aod | aie 


lq lq| laq|—1 


E\-£); Hy +2E, +(E, +£)); TEs, for s=n—2 
and under this equivalence V becomes O (mod 1) except for e = 
have V = (0, 1, 1)/2 (mod 1); 
or 


a 


E, — E, we also 


(ii) f is equivalent to a positive multiple of 


° ° 9 
ae eens ae 2 a 
Xe tt TX ab EX, PE or se 
ee 
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for suitable choice of ¢,, fg, and under this equivalence V becomes (1/2, v) (mod 1) for 


h and V given in Table V. 








TABLE V 
k h V (mod 1) 
0 + 2(E; + £1) 1/2 (ui, u2); 4, 56 uz (mod 2) 
+ 2(E; — £)) 1/2 (uy, 42); wu) = ue (mod 2) 
+4E£, + 1/4 
+2E, aj/4 (i 1; 1) 
4H, 1/2 (u;, U2); uy Ue = O (mod 2) 
2E4 1/2 (u + », u, v, 0); (u, v) FE (0, 0) 
+2E; + 1/4 (1, 1, 2, 0, 2) 
ae +4E, 0 
+ 2E3 1/2 (1, 1, 1) 
+2Es 1/2 (1, 1; 0, 0, 0) 
+2 -2(E, + £)) (0, 0) 
+ 2(E, + £)) 1/2 (1, 1) 
+ 2(E1 — E,) 1/2 (0, 1) 
A 
F2Es 0 
A 
+2Eg 1/2 (0, 0, 1, 1, 0, 0) 


ee 


Moreover, if P; (f)"/ | D| < 1/4, then PP (f)9) | D | < Ys. 


For the nonnegative spectra we have 

Theorem 3—If f is an indefinite quadratic form of rank n > 21, then Py (f)"/ 
| D{ «15. 

Because of the recent work of Margulis, the restriction 1 > 21 may bea tech- 
nical one, which is ‘used only to obtainthe bounds in the proof of Lemma 6. In 
particular, the work of Bambah et al., shows that the constant given in Theorem 1 
also holds for n > 2,s = 0, 1,2; forn > 7,s = —1. 

2. A CONJECTURE OF OPPENHEIM (ALSO CALLED DAVENPORT’S CONJECTURE) 

In Watson® the author defines the function G (f, V), for fixed V € R", to be the 
infimum of y for which the inequality 

cx fR¥+V set 7 
is solvable in integral X € R’, for allc. Hence, for all V € R", 


P(f,V) < P*(£,V) SGU, Y). 


Setting 
G(f) = sup {G(f, V): V € R". 
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Theorem 2 in Watson*®? shows that, for n > 21, G(f)"/|D| <1. each 
proves the same inequality for zero forms of arbitrary rank. More recently, at 
etal.®, proved that, for «4, B>0O, «+ 68=2|DJ|'", the strict inequality 
—a < f(X + V) < @ holds for any nonzero form f with | s(f) | <& 2. 
Restricting our attention to indefinite forms, we define 
M (f) = inf { | f(X) | #0: X is an integral point}. 


In Oppenheim” it was conjectured that, if M(f) #0 andn => 5, then fis a 
rational form. In Theorem 1 of Margulies*! it is proved that if fis an indefinite 
quadratic form of rank n > 3 which is not a multiple of an integral form, M (f) = 0. 
We shall see below that this result allows us to consider only rational forms. 


Lemma |—Let f be an indefinite quadratic form of rank n > 3. If fis nota 
multiple of a rational form, then G (f) = 0. Moreover, if fis rational form of rank 
n » Sand V is not a rational point, then G (f, V) = 0, 


Proor : We first suppose that fis not a rational form. Then, by Oppenheim’s 
conjecture, M(f) = 0. In Oppenheim?* it is shown that, forn > 3, P+ (f,0) = 0 
implies P* (—/, 0)=0; hence, from M (f) = 0 we obtain P+ (f,0) =0. In Theorem 1 
of Watson” it is proved that G(f) = 0 follows from P+ (f, 0) = 0. This proves 
the first assertion. 


If fis an indefinite integral quadratic form with n> 5, then by Meyer’s 
Theorem [Cassels’, p. 75, Corollary 1], f represents zero nontrivially. Hence, by 
Theorem 2 in Watson®, G (f, V) = 0 for any nonrational point V. 


3. RESULTS ON EQUIVALENCE 

In this section fandg denote nonsingular primitive integral (not necessarily 
indefinite) quadratic forms of respective ranks n (f), n (g), and respective signatures 
5(f),5(g). Because in any of our spectra the values for two equivalent forms are 
equal, it suffices to consider a complete set of representatives for the equivalence 
classes. It has been proved, for instance Theorem 1. ] in Cassels’, that the number of 
equivalence classes of forms of fixed rank and discriminant is finite. Because results 
to consider only a finite set of discrimi- 


Definition—For prime P, fis saidto be p-adically equivalent to & (denoted by 


tp g) if, for each positive integer ¢, there exists an integral matrix N such that p does 
not divide det N and I(NX) =e (X) (mod p‘). 
Definition—Let d (f) = d(g) = d. Then f is said 


to be congruentially equi- 
valent to g (denoted by f= g)iff pg, 


for all primes p dividing d. 
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Femma 2—Let f and g be indefinite forms of the same rank n > 3, and the same 
discriminant d. If, for any integer m > 5, dis not divisible by m""")/?, then f ~ g 
if and only if f= g and s (Ff) =s (g). 


PRoor: Watson®!, Corollary 2, p. 111. 


We shall later prove that all of our discriminants are small, relative to the 
hypothesis of Lemma 2. 


Lemma 3—Let pbe any prime dividing the integerd. Then there are only 
finitely many p-adic equivalence classes of fixed rank and discriminant d. Each class 
1s represented by a disjoint sum 


fo + Ph, +...+ ptf (4) 


where, for each O Ci <k, either /; = 0 or there exists t; > 0 such that; = A; +g, 
where either g = 0 or, if p = 2, 


g = E., + cx’, x? + cy’, —x? — cy* fore = 1, —3 


if p is odd, 
g = x®, x2 — by?, bx? 
where 5 is quadratic nonresidue (mod p). 
Proor : Watson®', Theorems 32 and 35, p. 54 fff. 


We note that since the sum in (4) is disjoint, the determinant of the form in (4) 
is det fo . det (pf,)..... det (p* f;.). 


Lemma 4—Let n(f) = n(g) and pbeany prime. Thenfpg if and only if 
det f - det g is a p-adic square, p‘|| d (f), d(g), and there exists an integral matrix N 
such that p does not divide det N and f(N X) = g (X) (mod p’). 


Proor: We note that the quotient det f - det g/(—1)*("/7] d(f) - d(g) equals 1 
or 4. This lemma thus follows from the analogues for determinants given in Watson®! 
the sufficiency is given in Theorem 33 (ii); the necessity can be found on p. 50. 


The last result will be used to eliminate replications in the list obtained from 
Lemma 3. To obtain the congruential classes for a fixed rank and discriminant d, the 
information from the p-adic classes for all p dividing d must be combined. 


4. A REDUCTION OF OUR PROBLEM 


We recall that in Lemma | we showed if the rank of / is at least 5, and either f 
is not a multiple of a rational form or V is not a rational point, then G(f, V) = 0. 
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We now return to the function G (/, V) and restrict our attention to primitive integral 
forms / and rational points V. 

Writing V = q7? (ui, «--, Un), Where gcd (u,,..., U,) = 1, expansion of f(X + v) 
yields f(X + V) = f(X) + I/c: L (X) + r, where r is a rational number and c1 is 
the integer such that Z is a linear integral form in which any common divisor of its 
coefficients is relatively prime to ¢c,. 


Lemma 5—gq divides c, D. 


Proor : For f (X) = >» a; Xi X;, expansion of f (X + V) yields the linear 
1<i<j<n 
form 
C1 
co 2 SS ous 
q 4g Ut Xy 
j—l im) 
where 


{° a, for i = j 
bi; a ; oa 
aj for ij. 


Setting N; = 5 bj, u, the fact that Z (X) is an integral form implies that 
Lid 


N; is an integer, for alll <j <7. Thus, forN = (N,,.... Na), U = (uj, ..., u,), and 
B = (b1,), we have gN = c, U B. 


Since det B = det f 4 0, then by Cramer’s Rule 


qB, 
Cl det f 


i = 





where B; is the correct minor of B augmented by N. We recall that U is an integral 
point with gcd (u,,..., u,, g) = 1 and so q divides c, det /. 


The following lemma follows from Theorem 3 of Watson’. Watson’s proof uses 
the restriction n > 21; in light of Margulis’ work a modification of Watson’s proof 
should allow a weakening of n > 21. As noted in Watson*, p. 567 the result may be 
true for all 7 > 9, but there is a counterexample for all n < 8; namely, 


f= xt 2 (xy + 1/2)? 4. + 2 (xn + 1/2)? 
has | det 4 | = 22"-1 £(X) = (s + 1)/2 (mod 4), for all integral X, 


Lemma 6—Letn > 21. IfG (fo V)9) iD |. > YS, thence 4 and one of the 
following holds : 


(i) G(f,V) = 1and| D| <5; 
GD)°G.Cf, V) = 251 Di. <— si 2", and f = L* (mod 2). 
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Recalling that G(f, V) > P* (f, V), > P (Jf, V), our analysis of each spectrum 
above 1/5 will be divided into the cases given in Lemma 6. 
5. PROOF OF THE THEOREMS FOR THE CASE G (f, V) = 1. 


We consider primitive integral forms f with n > 21 and rational points V = 
for which G (f, V) = 


U 


Q le 


Lemma 7—If n > 5 and P+ (f, V)"/ | D | > 1/5, then P* (f, V) = 1. 


Proor : Since G(f,V) > P*(f,V), then G (f,V)"/ | D | > 1/5 and the hypotheses 
of Lemma 6 hold. Thus, c, = | and from Lemma 5 we obtain that q divides D. Hence, 
the denominator of the rational number P* (/, V) divides D*. Hence, either 
Pog, Vy al or Ps (7,.V)'<. 1 —-1/D*. 


eee oe len FF (7, V) <.,1 — 1fD*. Setting J ( |D),n) = 
(1 — D-*)"/ | D| , we obtain P*(f, V)")|DI<T(|D|,n 


Sincen > 2l1and!1<|D|<4,7T(|DJ|,n)is a decreasing function of ” and an 
increasing function of | D|. Therefore, T(| D|,”) < T (4, 21) < 1/5, contrary to 
hypothesis. 


Lemma 8—Let t > 1 be an integer, p be an odd prime, and » be any quadratic 
nonresidue (mod p). If d, and d, are integers which are not divisible by p, then 
dy x3 + pid, xi pe (kx? + p' x4) 
where cd, = | (mod P‘) and k = 1 or b is chosen so that kd, d; is a p-adic square. 
Proor : For convenience, we define the forms 
g, (X) = d, x? + p' d, x2 and g,(X) =c (kx? + p! x3) 


for any c, k as given in the conclusion. Since kd, d, a p-adic square, there exists m 
such that kd, d, = m’ (mod p'). Also, 


det g, - det g, = 16 c® kd, d, p2! is a p-adic square and pl d (g,), d (g,). For 
N= oe, oy a where det N is not divisible by p, 
go (N X) = c(k (mk x,)? + pt (de x.)”) 
= cm? k™ x? + p! cd? x3, 


where 


cd, d, = d, (mod p’') 


Ill 


c m= k-* 
and 
cd? = d. (mod p’). 
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g. (N X) = g, (X) (mod p‘) and, by Lemma 4, g; p gp. 


Lemma 9—Let f be a primitive integral quadratic form withn(f) =n > 4. If 


| D| <4, then, for some gas givenin Table VI, f= H,, + g, where the sum is 
disjoint and 2m = n — n(g). 


TABLE VI 





| D| 1 2 3 4 





g 0 +5; +E, E,\ +E); £,-£,; 2H Es: +2E); +E> 
Moreover, these forms are pairwise incongruent. 


Proor: If | D| = 1, then the integralness of dimplies that nis even. Since 
n is even, then d = 0, | (mod 4) by (3) and sod = 1. Moreover, d (Hm) = 1 implies 
that f = Hm by definition of congruential equivalence. 


If | D| = 2andn were even, then | D| = | d| = 2 (mod 4), a contradiction 
tod =0,1(mod4). Hence,n is odd and|d| = 1. Also, d(H, + £)=+1. 
Therefore, choosing the sign so that d = d (H,, + E,), we have that f = H,, + E,. 


If| D| = 3, the integralness of d again implies that n is even andsod = — 3 
= d(+ E,) =d(H,, + E.). Hence, it is sufficient to show thatf~ Hm + E,. By 
3 


Lemma Fi Bes fo + 3f,. Also, Lemma 4 implies 3 || d( fy + 3f,) and son (f/f) = 1. 
Since n is even, n (f,) must be odd and 
J Hee Pid, XO ERs ee 
3 


for some d,, d, $ O (mod 3), Applying Lemma 8 we obtain 
Pe te Oe eee 3) 
3 


for cd, = 1 (mod 3) and k = + | chosen so that kd, 4 


? 2 1S a padicsquare. For the 
choice k = —- |, we have 


d(f)+ d(Hm + ¢(-x*_ + 3x?) = — 36e? 
which is not a 3-adic square. Therefore, 
Py He be Reese): 


11 
For N =| - ; rs we obtain E, (NX) = x? + 3x2. Since det N = 


=— 2, by Lemma 4 
I~ Hm + Ex, proving f = H,, + Ey, 
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For the case | D| = 4,f ~ fy + 2f:, Lemma 4 implies that 4 || det (/o + 2/1). 
2 
Hence, n (f,) < 2 and each of f, f; is of the form Hm + g, where 
g=0, E., + cE), E, + cE,, or — FE; — cE,, force = 1, —3. 


Moreover, 4||2" (/,) det f, - det/| restricts f, to0, Mi, FE., + cE;. We consider 
each of these cases in Table VII, where we let R = 0, 1, 2 be such that 24|| det /o. 


TABLE VII 
fy R fo det f + det (fy + 2f1) c fo + 2f 
0 2. A mictirl etCka +16c 1 Hm + E, + Ey 
Hm — E, — cE, +16¢ 1 Am — E, — E; 
A, 0 Hy, +16 — Ay + 2H, 
Hy + E, +48 — — 
E; 0 H,, +48 _ _— 
Hy + E: +144 — Hm + Es + 2E, 
+cE, 0 Hm ani Wey? 1 Hye 25; 
Hy + Es +48c —3 Hm + E, + 6E£, 





We next apply Lemma 4 to show that FE, + 2E, = E, and E, + 6E, = Es. 
We recall that 2'‘|| d (/), set g = fo + 2/;, and note that 





TABLE VIII 
g d(g) t N det N g (NX) (mod 2%) 
feed 0 1 
11-1 -1 
36 #4 —1 Ex (X%) 
gtts 071 oO 0 
ea 0 0 
jE Sit 
Eo + 6Ey 18 1 he ph a) 1 Es (X) 
O20 SI 


reer ——————E—E—E—EeEEEeEe——_—_—__ 


Hence, E, + 2E, = E, and £,+6E, = Es and Table VI follows from Table VII. 


Finally, we note that s (Hm + g) = s(g). Since the signatures of congruent 
forms are in the same congruence class (mod 8), to complete the proof of Lemma 9 it 
suffices to check that there exists no matrix N of odd determinant such that 2H; (N X) 
= (E, — E,) (X) (mod 4). 


Lemma \0—Let f be a primitive integral quadratic form withn (f)=n2>4 
and s(/) = k (mod 8), where -3 < k<4. If|D| <4, thenf = H,, + g, where 
Im =n — n(g), s(g) = k, and g is given in Table IX 
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TABLE IX 
a a i a 
k 0 +1 2 73 - 
g 0; E, — Ey; 2H; +E); +25, +Ex; + (E, + £1) +E; Ee 


Moreover, the forms are incongruent. 


ProoF: From Lemma 9 we obtain f= Hm +g, for2m =n — n(g) andg 
listed in Table VI. Since s (H,,) = Oand f = H,, + g, thens (f/f) = 5 (Hm + g)=s (g) 


(mod 8). For all g listed in Table VI we have —- 3 << k < 4. Hence, s(g) = k and 
Lemma 10 is proved. 


Proposition 1—Let f be a primitive integral quadratic form with n(f) =n => 5, 
|D| <4, anid s (f) = 8q + k, where -3 <k < 4. Asin Lemma 10, f = Hm + g, 
where 2m = n — n(g) is given in Table IX. Let 2¢ = n —n(g) — 81q]. 


(i) Ift > 0, then f ~ sgn (q) - Fig: + Hi + g, where g is given in Table IX. 
(ii) If t < 0, then f~ sgn (q) - Figi-it +h, where h is given in Table X. 


TABLE X 
a ee ee 
k t h 
A 
+2 ah H, + £6; Hi - Le 
+3 orl tae Aits3 + Es 
4 5h 5 Boa 


$e 


Moreover, these forms are pairwise inequivalent. 


Proor : Since equivalent forms are congruent and have the same Signature, the 
given forms are pairwise inequivalent. 


Applying Lemma 10 we have that f = Aaigias + 8, 


We note that since d(+ E,) = d(-+- Hi) = 1 andn 
definition we have + EF, = 


where g is givenin Table IX. 
(+ Es) =n (+ M,) = 8, by 
+ H,andsoforallr >0,4+7= + Hy. 


Considering the case when t > 0, we have J = sgn (9). Figt + Hy + g with 
§ (sgn (9) - Fit + He + 8) = 8q + s(g) = 89 +k =s(f). From Lemma 2 we 
thus have f ~ sgn (q) - Fig + A, + g. 


We next note thatis(f)|=8jq|=+ Ps(9) 1 IEF <:0e then q ~ 0 and 
RS Sse Cis (7) lee ls (g)|) 
n—s(f)— (2(g) +1 s(g)|)> 2—(n(g) + | s(g)|) 

Since fis indefinite. In particular, n (g) + |s(g)| > 4. Considering the choices of 
8, we thus have that | s(g) | > 2 and Is(f)|=8 lq|— Is(g)|. An inspection 
of the forms in Table 1X, combined with 21 — n(f) — n(g)—(l 5 

sae Cl s(f)1+ | s(g)]), 


-2 2 2t=n—n(g)—8\q| = 
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TABLE XI 
k g n(f)- Is(M|< d(f) t 
2 +E; + (E; + £)) 2 —3; —4 ced | 
+3 +E3 4 —2 —1; —2 
$ Es 6 4 —1; —2; -3 





Hence, when ¢ < 0,f = Ay gay: + g = sgn (q)- Figti-1 + Ar4a + 8, for g given in 
Table XI. 


To prove that f ~ sgn (q). Fi,t-1 + A, for the forms A given in Table X it thus 
suffices to show that 


Aisg + g Shands (h) =k — 8 - sgn (k) a4 5) 
the latter since we wants(f) =8 - sgn(g) -(|¢|—' +5 (A); 5 (4) = (8¢ + &) 
—8q +8 - sgn(q) =k + 8.sgn(q). We recallthat |8¢+k|=|s(f)|=8: 
|q| —|k|andsog <0 ifand only ifk > 0. Hence,k +8 - sen(q) =k —8- 


sgn (k), which gives the second condition of (5). 


We now consider the positive definite forms Eg, E,, E;, E, with respective dis- 
criminants —3, —4, 2,4. By Lemma 10, each of these is congruent to Hm, + g, for 
appropriate m > 0 and some g in Table IX. An analysis of discriminants, ranks, and 
signature (mod 8) yields 


A 
E, = H, — E,; Es = H, —(E: + £,); Es = Wi — Es; Ey = + Ey. 


The conclussons of the proposition thus follow from our above analysis and the infor- 
mation in Table XI. 


Proposition 2—Let f be a primitive integral quadratic form of rank n > 21 and 
let V € R® for which G(f, V)"/ | D| > 1/5. If G6(4 VY) = 1, thnf~ + F,+ H, 
spre for the appropriate choice of sign, and r, Y, a given in Proposition 1. Under 
this equivalence V becomes (0, v) for Vv © Rs), Moreover, {f(X + V): Xis an 
integral point} = 2 (V) + Z. 


Proor: Since G(f, V)"/ | D| > 1/5 and G (f, V) = 1, by Lemma 6 we have 
| D | <4 and Proposition | implies that f~ -— Fr + H, + g, for the appropriate 


choice of sign, andr, v,g. We write V = q™) (ui,..+5 Un), Where gcd (u;, ..» Un) = 1, 
and as in Lemma 5 f(X + V) =f (X) + L (X)/c, + r, where Z is an integral linear 


form and, by Lemma 6, c; = 1. Since + F, +H, + isa disjoint sum and the 
statement of this proposition only gives information about + F, + A,, it is likewise 
sufficient to consider only FE, and H,. But 
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Hi (x + ulg,y + u']g) = M(x, y+ qx + uy) + uu'}q? 


and the linear portion is an integral form; for all 8r < i < 2m, v) = O (mod 1). 
Expanding E, (X + V) we obtain that the linear portion equals 


1 4 4 is 
q | = (2a + P (uy — Bi) 1 + tos (2144 Pr u,) xive- | 


Since this is an integral form, 


q divides 2u; + 2 (uj — Uj4,) ...(6) 


and 
q divides 2uj,, — 2 uy, ...(7) 
#4 


for all i. Summing over i in (7) we obtain 
q divides 2U, — 3U, ...(8) 


where 


4 4 
U,= © uy, andU, = 2 wy. 
i=] f—1 
We recall from Lemma S that q divides ce, D= Dand|DJ|< 4. Therefore, 
2 and 3 are the only prime divisors of g. 
If 3 divides g, then (8) implies that 3 divides U,. Adding (6) and (7) 
q divides 2u; — U, + 3ui,4 


and so 3 divides each u;. Hence, from (6) we obtain that 3 divides each U, — wi44. 
Therefore, if 3 divides qg, then 


uj; = O (mod 3) for alli < | 8r |. 
If q is divisible by a power of 2, then (8) implies ~ that U, is even. Subtracting 
(7) from (6) we thus have that U, + uiy4 is even; 


4 
x (U; + w44) = 5U;, is even; 


and so w,,4 is even, for alll <i <4 The evenness of ui, for i & 4, follows from (6) 
and the evenness of each wj44 ae U;. Hence, there exist integers w, such that 
uj = 2), for alli < 8. If 4 divides q, then (6), (7), (8) can be rewritten in terms of 
the w,; the above argument applies and yields w;, = O (mod 2), for all i < 8. 


Pe Therefore, u; = O(mod q) for all i< 2m and we write V = (0, Vv) where 
VE RX), 
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To complete the proof, we show that 
(f(X + V): X is integral} = g(V) + Z. ...(9) 


Writing X = (Y, Z), where Z € Zre), we have X + V=(Y,Z + v). We recall, 
from Proposition 1, that 7 (g) S 6 and so| 8r|+ 2» > 15, implying that either 
|r| >Oorv> 0. 

1.020 then {(+ F, +.) (Y): integral Y}:= Z and, for fixed Z, {f(Y, Z + v): 
Y is integral} = a (Z + Vv) + Z. G(f, V) = 1 thus implies (9). 

On the other hand, if v = Oand s (+F,) > 0, then {+ F, (X) : integral X} = 


the set of nonegative integers. 


Therefore, for fixed Z € Z2), {f (Y,Z + V): integral Y} = g (2+ V)+4+ N*. 
From G (F, V) = | we thus obtain 


&§(Z+ V) — g(V) E Z, for all Z € Z"2), 


Moreover, since f is indefinite and {F, (Y): integral Y} = IN*, there exist 
integral vectors Z; such that g (Z; + V)—+ — cc. Writing g 2 (Zi Vy Sav 


i>oo 
= m, we have, for fixed /, {f (Y,Z; + V): integral Y} = g 2 (V) +m + IN*, from which 
(9) follows. An analogus argument applies for s (-- F;) < 0. 


Proof of the Theorems of the case G(f, V) = 1 
We assume ) that cae Oa Pa i | = | > 1/5 and G(f, V) = 1. By Proposition 1, 
f[~Fr+H, + g Pe foe some r, v and . as given in Esopoattion 1. Moreover, from 


Proposition 2. under this equivalence we have V = (0 v). We claim that all non- 
integral V which satisfy G ( f, V) = 1 are given in Table XII. 


Column 4 (Table XI!) follows from Lemmas 5 and 6. We shall complete the argu- 
ment fare = + E;. Recalling (9), we obtain that, for any integral N = (m, Mm, 73), 
E, (V + N) — EB; (Vv) is an integer. Setting V = 1/4 (u,, u., us) we thus obtain that 

u, (2n, + ny, + 3) + uy (m + 2n, +73) + Us (Ny + My +2n3) = O (mod 4) 
for all choices of integral N. Using N = (1, 0, 0), (0, 1, 9), (0, 0, 1) we thus obtain 
that uv, = iw = uz = | (mod 4), completing the argument forg = + Es. 


Continuing with the proof of the theorems for the case when G (f, V) = 1, we 
assume that P+ (f, V)"/D | > 1/5, where 
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TABLE XII 
k a | D | q#l V (mod 1) 
0 1 i — 
1 
E, — E,; 2H, 4 2,4 zy (4, ue) 
a 
+1 Ey 2 2 2 
u 
2Ey 4 2,4 +f 
u 
2 E; 3 3 3 (1,1) 
1 
Ey “2 Ey 4 an ot a (ui, U2) 
u 
Es 3 3 3 (Onis 1. Fo Osean) 
A 1 
Es 4 2,4 > (0, 0, u, v, u+yv, uty) 
u 
a8 E3 a 2,4 4 dG, 1, 1) 
u 
Es 4 2,4 iz Cie 22 9) 
1 
4 Ey 4 2,4 > (ut+y, u, v, 0) 
ec ee ere 
TABLE XIII 
fra: ~~ ~~ 
k Te tg (V) V (mod 1) PT (f, V)"/|D| 
0 _- 0 1 
i re 2 u 
Ei, — £, 4 4, —u ) 2 I, 1) 1/4 
Ay 4 (u1 us) 0 1/4 
zee Ey u?!4 0. 1/2 
2E u2/8 0 1/4 
FE, + E, + (ud + u?) 0 1/4 
E¢ 2 u2!/3 0 
- 1/3 
E6 £ (u’+y2) 0 1/4 
= ah | E 9 
3 3 u7/8 0 1/4 
E 2 
5 5 u*'8 0 1/4 


i) 


af Eq 4 (u?+v? + yy) 1/4 
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P* (f, V) = inf {fX + V) > 0: integral X}. 
By Lemma 7, P‘ (f, V)"/ | D | > 1/5 implies that 


Pt+(f, V) = 138 (V) € Z, by (9). Using this and Table XII, we obtain column 5 
Table XII{ below and so Theorems | and 2 for the case when G(/, V) = 1. 


We note that P(f, V)” | D| > 1/5 implies that Pt (f, V)"/| D| > 1/5. There- 
fore, ifG (f, V) = 1 and P (f, V)"/ | D | > 1/5, then the above argument applies with 


V and g as listed in Table XI. Since 
P (f, 0) = inf { f(X) > 0: X is nonzero integral point} = 0 


for any of these /, then P; (f)" | D| < 1/5 for all fand Theorem 3 is proved for the 
case when G (f, V) = 1. 


6. PROOF OF THE THEOREMS FOR THE CASE G (f, V) = 2 


By Lemma 6, we have that f = L* (mod 2) and so, fori + j, aj; is even. By 
definition, the diagonal entries of the matrix A are even. 


Lemma nai g =fy + 2fi, then fo = cx? or + x? + cx} . 
: 2 


Proor Since f 3 g, by definition there exists an integral matrix M with odd deter- 


minant such that g (M X) =f (X) (mod 2). A comparison of coefficients yields that, 
for yi = Mi; X1 +... +Min Xn, fo (Y) = fo (X) (mod 2). If fp = E, or there exists t > 0 
such that /, = A; + go, then fori #/ 


My My + M,, My = aj; (mod 2). 
Since a; is even, then also 
My M2; — M,; my = O (mod 2). 

Letting M,; be the n—2 order minor of M with the first two rows and the j, j 

columns removed, then (cf. Hadley’’, page 98, eqn. (3-100)). 
det M = Se, (my m2y — m,; my) det Mi;. 
for the correct choice aha jy = +1, contrary to the oddness of det M. 

Lemma 12—Let n> 5 and let aj, be even, for alliA~/j. If i | See) en 
f~ + xi +g, where the sign can be chosen so that g is an indefinite form of rank 
n— 1. 

Proor : We first show that the result follows if f represents both + |. Let B be 


an integral point for which f(B) = + 1. Then Bis primitive and B can be extended 
to a basis, thus obtaining a unimodular matrix B in which B is the first column, Thus, 
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f (BX) = 4 xt 4X1 M (Xe, gg Xp) 1 A(X g; cee Xu), 


where A is a quadratic form and M is the linear form given by 
M (&) =) So2at BOS be a on Be ee ee 
i=] k=2 i<j k=2 
Since a; is even for all i + j, then M has even coefficients, say 


M = 2 (cg xt + ... + cy Xs). 
Hence, 


S(BX) = + (%) E(x. +... +c, x,))? + (h F M?/4) (x2, ..., Xn)? 


f~ + x) + g, for the quadratic form g = h = M?/*, The sign of + 1 is chosen so 
that g is indefinite. 


In order to complete the proof of the lemma it suffices to show that f represents 
both + 1. For this weshall use Theorem 51 (ii) in Watson®*: Let a be a fixed 
integer. If there exists a primitive solution X to f (X) = a(mod d), then there exists 
a form h which is congruentially equivalent to f, has the same signature, and properly 
represents a. 


We note that | D | < 2"*? and Lemma 2 imply that such / is equivalent to f 
and hence that f itself properly represents a. Moreover, from jf = L* (mod 2) we 
obtain that 2" divides D and so | D| = r. 2", for somer < 4. 


For the case when r = 3, by Lemma 3 we write tf ~ go + 3g,. Lemma 4 im- 
3 


plies that n (g,) = } andson (g,) > 4. Inspecting the choices in Lemma 3, we observe 
that g, = H, + g, for some t > |. Therefore, go represents + |; that is, each of f (X) 
= + | (mod 3) is solvable. Suppose we have shown that, for each of « = + 1, there 
exists primitive P such that f(P) = « + 2"m. For primitive Q such that 


f(Q) =« (mod 3), weset P =P + Q. We note that, for the associated bilinear form 
defined by 


4 (X, Y) re Qj Xt Vj, 
f (P,P) = f (P) + £(P, Q); 


J (P) =f (P) +2¢(P, Q) +6 (mod 3). 
Therefore, for T = 2” 2 (mod 3), 
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f (P+ 72-1 P) = f(P) +2" TLR, P) 4 22 T? f (P) 

=f (P) + 2"(T (fF (P) + /(P, Q)) + 2" (f (P) + 2F (P, Q) 
+ €)] (mod 3:2”) 
= f(P) + 2° [2°72 f(P) + 2"-? ¢] 
= € + 2" m+ 2"[2"1 (e + 2m) + 2%] 
Sse + (2" + 2°-*) m 
= ¢e (mod 3:2"). 
Hence, for any | D| < 2"*° it suffices to shows that each of 
f (X) = + 1 (mod 21) is solvable for 2'\Id. 


By Lemma 3, we write / wee: -+ 2g, = hand consider the choices fort. al, 


for the some v > 1, f= A, + fy, then fy properly represents both + 1. We may 
thus assume v = 0. Also, since fis primitive, /) represents some odd integer. There- 
fore, if g1 = Hy, +8, for some w > 1, we would obtain that A represents all odd in- 
tegers. We may thus assume that v = w = 0. We write g, = fi + 2g. and observe 
thatn (fo),n (fi) <2. From Lemma 11, /; a E, and so 2"(fo) divides det f), and 
2"tne) divides D. Since | D | <2"+*, we obtain that g. = H,, E., exe anes, | 
= 2"? f, = E,. Also, since h ~ f, then det h' det f is a 2-adic square; that is, its 
2 


odd part is congruent to | (mod 8). Letting d, be the odd part of det g, and recalling 


that fo = bx} or + x? + bx} we thus obtain that 


bd, = + 3 (mod 8). ...(10) 
We also note that, for g = — xj + 2 E, (x, Xs). 
n (1,0) = — land g(-—1, 1,0) = 1. Hence, we may assume that b+ — 1. In 
Table XIV below we obtain Pi, P, such that A (P,) =1 (mod d) and A(P,) = -— 1 


(mod d) for the remaining choices of h. The congruence in (10) andn > 5 are used 
to obtain the possibilities for fg given in column 3. 


Proposition 3—Let n > 4. If fis an indefinite primitive integral quadratic form 
with for each i ¥ j, a); is even, and | D| <2"*?, then there exists r, such that f'~ + xj 
+... x* + 4h, for suitable choice of signs, where sum is disjoint and / is given in 


Table XV 
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TABLE XIV 
Ss SSS 
82 d, fo n | d | P Po 
E, 3 xe oi poner (1, 0) (1, 2,1, 2, 2) 
= + x 6 2° di, 0) di, 0, a; LB “3 2) 
A, l 32. 5 25 (1105+ 18) (1, 0; 0, —1, 4) 
— 3x, 5 26 (1, 0, 0, 1, 1) (1, 1, 0 
x 3%... bana (1, 0) (0, 1. 0, 0, —1, 1) 
op Ek ee (1, 0) (0, 1, 1, 0) 
=e I Ce ck Oa 28 (1, 0) (2, 1, 2, 4, 0) 
£;=Bagy ces (1, 0) (0, 1, 1, 0) 
3x" 3 = au Ae 6 
: ee | 2 (1, 0) (1: 0. "Te 44) 
= 3x" 3 “s +. 5 6 
x. x, x, 5 2 Ch. 0) (1, ay 1) 
LL neers 


TABLE XV 


1 0 
2 2Ey 
3 2E.; 3E 5 


| PROOF : Since, for each i ¥~ Jj, ai; is even, then Qn divides Dand | D | < 2+ 
implies that | D | = r°2", for some | & a4: 


Ifn > 5, we can use Lemma 12 repeatedly until we obtain 
a EE ee ie sa 


for some indefinite integral quadratic form Ji. The process stops when either n(f,) = 4 
orf: is not primitive. We set ny = n(f;) and D, = det f,. 


We first consider the case when f; is not primitive and n, = 4.Since | D, | = 2%r 
then f, = 2f, for some primitive indefinite quadratic from f, with n (J =m>4 
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and | det 2 | = r. Recalling that Lemma 10 and Prop. 1 were obtained without using 
the condition G (f, V) = 1, we thus have that f, ~ + F, + H, + g, for some choice 
of sign, vy, t > 0, and g given in Table IX or X. Recalling that + F, = Hi, and 


(—x,) (2) = — x1 x, and (x, + x, + Xa)? — 2 (x1 + x3) (x. + x2) 


2 
Xs 


a 


- ~ 


+x 


ts 09 


we obtain that 


+ xP +2(4 4) ~ xt ... k xX; 44» for some choice of sign...(11) 


Thus, for the casewhen /, is not primitive, each h = 2g appears in Table XV. 
Therefore, it suffices to consider the case when 
| ncn 2, Site ee a a SP oe 
where /, is primitive and | D, | = 16r. 


If r = 3, we first analyze 3-adic equivalence. From n, = 4, 3 ||Di, and Lemma 3, 
we have f, = H, + cx} + 3dxi. Since 


\ 


Hy, (%;—X2, X, + X2) = x} — xg and cx, + 3dx? = + x3 + 3x? (mod 3) 


then 
Gee mee Xe EX, PP Sxe 
3 
that is, there is exactly one 3-adic equivalence class for each of D; = 48, D, = — 48. 
D, = — 48. Hence, for r = 3 (and so for all |! <r < 4), congruential equivalence 


follows from 2-adic equivalence. We write fi 26 + 2g,; Lemma 11! implies that 
fe = CX, OF ge = x, + cx, . 


Ifr = 1, then det g,is odd andso g, = E, or g; = Aj. Since n, = 4, by 

: 9 Pa 2 2 44 : : 
benma 4 either s, = 2 My 3x +.2B, or fi Se tala) Hy, with the latter 
equivalent to + x? +... + x; by (II). We note that 3 (x1 + x2)?~(%1 + 2 + Xs) 


— (x1 + x3)? = x + 2E, (mod 4) and so Lemma 4 implies that 


ral ke) 
x, + 2E,~ Briere Xp nt ky : 
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Since + 3x) + 3x, = + x a ee (mod 4), again using 


Lemma 4 we have + 3x? + 3x, ~ + x; + x5. Combining these results, 


& xi tt 3x, + 22, ~ 
2 
saint te ee Ske et ike 2 4a! (oy ey a 
that is, = 0 whenr = 1. 


In a similar manner, if r = 3 then either /, ~ + xf +x} + 2E, 
2 


orf, ~ x; + 3x} + 2H,. Again using (11), (12), +3—= = 1 (mod 4), and 
2 

Lemma 4, an analysis of disjoint summands yields h = 2E, when r = 3. 
Recalling that 1 ~ go + 2g, with gy = ex! or + ie + CX, 

for r = 2 we have 2\/ det g;. Fromn(g,) > 2 we thus obtain g, = H, + bE,. 


Hence, /; = cx, + 2 (x2 x3; + bx} ). By considering —f, if necessary, we may 
assume thatc = 1,3. Also, since + chis a 2-adic square, we may aSsume that 


b= +c. Ifc =1, then (11) implies that fy ~ xy ty oe hk a 
For c = 3, we consider the cases b = + 3 separately : From 


(x, + 2x)? +2 (x1 + x,)? = 3 (x + Dice ) (mod 8) 


and Lemma 4 we obtain 3 (x; + 2x3) ~ x? 4 2x0 
2 
Hence, using (11), 


2 2 2 
3 (x + 2xq) + 2H ~ xt + 2x + 2H, ~ x8 + xt — x? +2x!; 


h = 2E,. Finally, for b = — 3, we observe that 
3 (x, + x2 + x3)? +2 (—%1 — x3) (x2 + x) 


= — Xt, Oe nd 4); 
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3E, + 2H, ~ — xi — x, — 3x} , by Lemma 4. 


Therefore, 
3x; + 2H, — 6x° ae Xi — xX) —3 (x2 + 2xt ) 
as above, andj = — 2), 
For r = 4, we again write /, ce &o + 281, where 2 = CX, or + xi + cx}; 4 


\|\det g,. Considering — f, if necessary, we may assume thatc = 1,3 and obtain the 
possibilities given in Table XVI, (using the determinant argument from Lemma 4). 
The empty entries of column 3 are discussed below. 





TABLE XVI 
80 &1 fi~ 
zy Hy 2x: xi ct x) — xy & 4x, , by (11) 
3x" H, + 6x, 
+ <: + “a + x, + 7 + x, + %, at 2x, se 2x! 
+x + 3x, + xe + 3x! 
+ aa iF cx, 2H C=) + i + cs + 4H, 
2EF,3 c= 3 


ER eer rere 


Each entry in column 3 of Table XVI is listed in Table XV. Noting that 
2 
3x2 + 12x} = 3x, F 4x} (mod 16) 
and 
3x} + 2H\= — x} + 2H, (mod 4) 


we obtain 


3x, + 2H, + 12x ~ 3x? + 2H, ¥ 4x; 


~ — xi + 2M F 4x4 


” 


2 2 2 3 
Mp Agee hg AX 
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the latter by (11). It thus suffices to consider 
lat Xb 3xy et exe oe 
and 
I ee xi + 3x, + 4E,. 
We note that 
cE (3x2)? + 3 (3x1 + 4x. + 2x5)? + 6 (x, + 2x2 + x3)? 


= xi +x, + 2x? (mod 8); 


2x1 + X_ + 2x3)? + 3 (x, + 2x,)? —6 (2x; + 3x2 + x3)? 


—xi — x, — 2x} (mod 8); 
— (%1 + 2x2 + 2x5)? + 3 (x, + 2xs)* — 6 (x1 + x, + X3)* 
=x} + x} + 2x% (mod 8) 
and also that the determinant of each transformation is odd. Hence, by Lemma 4. 
+ xy + 3x2 go Ox er + x, + 2x? x! 


for some choice of signs and So the fourth row of Table XVI hash = 2(E, + E,) 


Moreover, 
3 (x, + 2x2 — 2x,)? + 4E, (2x, + X_,— x3, x, + xs) 
= —xt+ 4x. x, (mod 16); | 
—3 (x, + 2x, + 2x9)? + 4E, (—x,, + Xs, x; + 2X2 + Xs) 


= xi + 4x, x3 (mod 16); 


where again each determinant is odd. By Lemma 4, + 3x? + 4F Pe x. ee 
’ 1 5 lie 1 1 


and the sixth row of Table XVI satisfies the conclusion with h = 4H, 
Finally we consider the seco 
nd ; ~ 
row; namely, /, : 3 (xi +4x,)? + 2A. We note 


that 3 ((x, + 4x4)2 +4 (x, + x,)?) = — x _ 4x, (mod 16); 
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3 (3x1 + 4x4)? — 4 (x1 + x4)') = — xi + 4x} (mod 16). 


Since each of these determinats is odd, 3 (x; + 4x,)> ~ — xs + 4x 


Hence, by (11), 


= 
i 


. 2 2 2 2 9 ° 
fey plies ee seek eX Lee 4X, 


1 


4E,. 


Summarizing, we are considering f for which f ~ -+ a ee eee tt efor 


some choice of signs, where /; is indefinite, n(f/,) = 4, and | det(/i)| = 7'2", for 
some r < 4. We have shown that there exists a choice of signs such that f, = + 


x? ,tx)_, +g, where g is given in Table XVII. 


TABLE XVII 





r 1 2 3 4 
g 0 2E, 2E, 4E,; 2(E, + E)); 4H; 





Setting h, = + x>_, + x°_, +g, we thus have that f; = h and sos(f,) =s (h,) 


n-3 na-—2 


(mod 8). Since /; is indefinite and” (fi) = 4, then | s(f1) | € 3 andsos(f,) = s (A). 
Lemma 2 thus implies that /, ~ /,, which completes the proof of Proposition 3. 


: [ies i. oe, o 2 3 
For convenience, we set Ki 510 = — xX) de 2x, = i x, +1 be ot AE Xt," 


Proposition 4—Let f be an indefinite, primitive, integral quadratic form of rank 
n>2l. If V € R" for which G (f, V) n/|D| > 5, and G(f, V) = 2, then 
f — Kiyst2 + h, for some h given in Table XV. Moreover, under this equivalence, for 


alll] <i<t, + t., % =4(mod 1). Also, 


(Kit, (X + V): X is integral} = Kr... (1/2) + 2 2. ae 3) 


Proor : From Lemma 6 we know that f = L® (mod 2) and so in particular, 
fori # j, a1, is even. Hence, we obtain the first conclusion from Proposition 3. By 
the disjointness of the sum it suffices to complete the proof for f = Ki,, For any /, 
we takeO < v, < 1. Again by disjointness, G(f,V) = 2 implies that (1 + vi)? 
— vi is zero or is at least 2: that is, % = 3. The set equality (13) follows from 


Fy 


Theorem 2 (i) in Watson**, completing the proof of Prop. 4. 
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2-—n|D| 
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TABLE XVIII 


PU pS a 
if P= (fF, V)= 2 
—t,*t.(mod 8) s(f) (mod 8) 


+h (V2) (mod 2) 








+h V2 (mod 1) 
0 ef 
u 
2E, al 
u 
2E, 3 Lhe 1) 
u 

2E6 3 (0, 1,1, 1,0, —1) 
2E:+E;) 4 (uy, us) 
2(E, - E;) 3 (uy, U2) 

u 
4E, 4 
u 

2E3 as tea 

4H, 4 (uy, uz) 

2E4 3 (u+yv, u, v, 0) 

u 
2E5 751, 2, 02) 
A 
2E¢ 


0 


0; u=0 (2) 
4; u=1 (2) 


tO 


33 u=--1 (3) 
0; u=0 (3) 


2 
mera. Ue u=+1 (3) 


0; u=0 (3) 
0; u), u2=0 (2) 
1; 4, W.=1 (2) 


4; 4. 56u, (2) 


0; u, =u, (2) 
— 4; u,=0, u,=1 (2) 
3; 4=1, u.=0 (2) 
0; u=0 (4) 
4; u=+1 (4) 
1; u=2 (4) 
0; u=0 (4) 
3; u=+1 (4) 
1; u=2 (4) 
0; uw; u.=0 (2) 
134; ue=1 (2) 
l;uAv(2) «+ 
0; u=v=0 (2) 
1; u=v=1 (2) 
0; u=0 (4) 
>, u=-+1 (4) 


1; u=2 (4) 


3 (0, 0, u, vy, u+y, u+v 0; u=v=0 (2) 


—4; uv (2) 
1; u=v=1 (2) 


oO 


—— 


HoH ot 
- N 


+ 
> 
N 
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TABLE XIX 
Pe ee hence. 
+2 (£,+E£)) 3 (U1, 42); u.Au, (2) 1/4 
+2 (E,-E)) 4 (ui, U2); u,) =uy (2) 
+4E; +} 
+2E, +11, 1,1) 
4A, $ (uy, u2); uy u,=0 (2) 
2E, 4 (1, u, v, 0); uv (2) 
4 (0, 1, 1, 0) 
+2E; +} (1, 1, 2, 0, 2) 
sl +2E) 0 1/2 
F2E, k 
+4E, 0 1/4 
F2E; 4 (1,1, 1) 
+2E; 4 (1, 1, 0, 0, 0) 
= a +2E; 0 1/3 
+2E.% 0 
+2 (£,+£)) 0 1/4 
+2 #2 (E,+£,) 4 (1, 1) 1/4 
+2 (E,—E;) + (0, 1) 
FIL 0 
+2 i (0, 0, 1, 1, 0, 0) 
+3 F4E, 4 1/4 
+2E; 0 
2k. 0 
4 4H, 4, 1) 1/4 
2E4 0 
2Eg 1 (0, 0, u, v, 1, 1); uv (2) 
TABLE XX 


Consider the case when P* (f, V) <2 





+h P*(f,W< Pt(f,V)n/ | DIS 
: a Ae 
23 3(3) 
2E,; 2E¢ 12 3 \ 24 
: 2(4)" 
Otherwise 4 I\38 


——— el 


958 MARY E. FLAHIVE 


Proof of the theorems for the case G(f, V) = 2 


If P+ (f, V)"/ | D| > t holds, then G(f, V)"/| D| > ¢. From G(f, V) = 2, 
we obtain f ~ Ki,,,. + A, for some choice of 4), t, and for some A listed in Table XV. 
Moreover, separating V = (Vj, V2) as in the sum Ki,,4, + h, by Prop. 4 we have Vi 
== (1/2) (mod 1). Also, from the disjointness of K;,,.. + 4, Lemma 6 implies that if 
M (X) is the linear part of h(X + V,) then / is integral and M* (X) = hA (X) (mod 2). 
Each / given in Table XV has even coefficients. Hence M?(X) = h X (mod 2) is 
equivalent to reguiring that each coefficient in M is even. Since A = 2g, for some zg 
ia Table XII, we thus obtain the first three columns of Table XVIII from Table XII. 


Column 4 is computed directly from the preceding columns. For the remainder of 
Table XVIII. we recall (13), obtaining 


L(V +X) = Kisty (1/2) + V2 (mod 2) = 3 (— tf) + te) + A(V2) 


+ h (V2) (mod 2). 


Hence, P* (f, V) = 2 exactly when — ¢, + 1, =—4h (V2) (mod 8). For the last 
column we use s (f/f ) = — t, + ts + s (+h) (mod 8), 


Summarizing, if P+ (f,V)"/| D| > 2, and Pt (f, V) = 2, then we have 
ft ~ Kis, + h where hand V = (0, V.) are given below (Tables XIX and XX). 


Since n > 21, each of the entries in the third column is less then ?. This com- 
pletes the proofs of Theorems | and 2 


For Theorem 3 we note that P if, V)"/| D| > 1/5 implies that P+ ChVAE D1 
> 1/S. Therefore, when P(AVY/|D|> 1/5 andG (f, V) = 2, the above argument 
applies and so P* (f, V) = 2. Hence, f and V are givenin Table XIX. But, for any of 
these we have P(f, V) = Oand so P(f, V)"| | D| < 1/5 for all f, V. 
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ON COMMON FIXED POINTS IN METRIC SPACES 


BARADA K. Ray 
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Some fixed point theorems for certain contractive type mapping are presented 
in this note. 


Throughout this paper (X, d) will denote a complete metric space unless other- 
wise stated and R*, the set of non-negative reals. Recently Kiventidis! proved the 
following : 


Theorem TK|—Let T be a self-mapping of X such that 
d (Tx, Ty) < d(x, y) — W (d(x, y)) ¥ x, y in X aL) 
where W : R* > R* is a continuous function such that 0 < W(r) < r for all re E€ Rt 
— {0} 
Then T has a unique fixed point : 
In what follows first we prove a theorem which gives Theorem TK | as a special case. 


Theorem 1—Let T be a continus mapping and 7,, T, be any other two Mappings 
of X into itself such that 


IT; = T;T Gi = 1, 2) a2) 
U* T(X) C T(x) 


and 


for all x, yin X (3) 


d(T, x, T: y) < d (Tx, Ty) — W (d (Tx, Ty)) 


Where W: R+ + Rt is a continuous function, with 


0< Wr) <r for allr € Rt — {0}. 


Then Fr,r,,7, = [YE X: x= Ty = Lik 


T, x} is non-empty. Furthermore 
Fr, spe Fr, = FT,T,,7; = 


{u}, for some wu in X. 


PROOF: Let Xo be an arbitrary point in Y. 
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Since 7T,(X) and T,(X) are subsets of Rix} we lete Tye, = Txs),i-and Te Aq 
= Txonyo,n = 0, 1, * Shae 


Then from (3) we have for alln > 1,x € pA 


= w (d (Tx,, Txr41)) < d (Txo, Tx,). 
r=0 


So the series of non-negative terms 


= Ww (d (Tx,, Tx,,1)) is convergent. 
[ sae 
From this it follows that lim W (d(Tx,, Txr4:)) = 9. 
roo 
Since W (0) = 0, so from the continuity of W we get 
lim W (d (Tx,, Txr,,)) = 0 


r+ 


> W ( lim d (Tx, TXr45)) = 0 


r->oo 


=> lim d (Tx;, Tx:4,) = 0 


roo 


which implies that {7x,} is Cauchy and so it converges to a point u in X, since X is 
complete. j 


PMereIOre 41 Xen 54) i Xanty (i Xencers= La Sands} And {71 Xo, = LeXen-1} being 
subsequences of {7Tx,} converge to ualso. But TT = T, T, i = 1, 2 and the contin- 
uity of T implies that Miss ffi 0 Pe) as ge 2 T(TXs,.;) Ta} lim” T, (7x) 


no n->oo 
= lim T (7:x.,) = Tu 4 fase Claas) ‘fe T (14X41) Tu 


n->oo n->oo no 


Now from (3) 
d(T, (Txon), To u) & d(T (Tx,), Tu) — W (d(T (Txon), Tu) 
Proceeding to the limit n —~ c°, we obtain Tu = Tyu. 
In a similar manner we can show that Tu = 7,u. Suppose u # Tu 


Now 
d (Ti (Tx) ToXon+41) 


< d(T (Tx0n), Txn41) — W (d (T(TXon))s TX2n41). 
Proceeding to the limit m — o°, we obtain 
d(Tu, u) < d (Tu, u) — W (d (Tu, u)), 


which is a contradiction. Thusu = Tu. 
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Sou = Tu=T,u = Tu, So Fr,r;,r, is nonempty. It follows easily from (3) 
Fr, = Fr, = Fr,r,,7, = {u}. 
Remarks : Theorem 7K | follows from Theorem | is one takes T 1 = T, and 
T = Ix where [yx is the identity mapping on X, 
In what follows we don’t take X as a complete metric space. 


Theorem 2—Let T be continuous mapping of a metric space X into itself satisfy- 
ing (1). If there exists a subset M of X and a point x, in M such that 


dl (x, Xe) —-d (Tx, Tx,) > Fd (x. Tx.) tor everyxinx — M ...(4) 
and if T maps M into a compact subset of ¥ then there exists a unique fixed point 


of T. 


Proor: Since T maps M into a compact set, Theorem 2 will follow from 
Theorem 7K | if it is shown that x, € M for every n, where x, = J" xo, n = 1,2,3.,... 
Let us suppose that x, + Tx 9. Then it follows easily that the sequence {C,}, where 
C, =Tdwoncten,. 1s non-increasing and since x, 4 TX), we get d (Xn, Xna1) 
at i (Xen Dk), 


But 

Pi ee oe < d(x,, Kati) Hh OL toe, 
So 

d (Xn, Xo) — d (Tx,, Txo) < d (x,, Xnai) + d (Xo, Tx) < 2d (xo, Tx). 
Hence it follows from (4) that x, € M for every 7, 


This completes the proof of Theorem 2. 
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Necessary and sufficient conditions for the asymptotic stability of a system of 
differential equations of dimension at most 4 with complex coefficients are 
established. When the coefficients are real, these conditions coincide with the 


well known Routh-Hurwitz conditions which we extend to include stability 
and uniform stability. 


l. INTRODUCTION 


Consider the homogeneous, first order linear system of ordinary differential equa- 
tions of n-dimensions X¥’ = A X where A is ann x nreal or complex constant matrix 
and X (t) is a column vector of the n dependent variables. The characteristic equation 


is a polynomial equation of degree n whose roots real or complex are the eigenvalues 
of A. 


We will discuss the stability of the system X¥’ = A Y when 4 is either a2 x 2 
3 x 3 or 4X4 real or complex matrix, and for this we follow the definitions of asym- 
ptotic stability, uniform stability and stability as given in Jordan and Smith*. 


First We note that when A is constant, then stability implies uniform stability 
(remarks following Definition 9.3 of Jordan and Smith’). 


Theorem 9.3 of Jordan and Smith® relates the question of stability to the nature 
of the eigenvalues of A. Details for the case of complex A can he found in Boyce and 
Diprima*. In case of asymptotic stability, we introduce the notions of Hurwitz poly- 
nomials and positive functions which will be defined in the next section and used to 
establish necessary and sufficient conditions for the asymptotic stability of the system 
X’ = AX where A is a complex matrix of dimension at most 4. When A is real, these 
conditions will automatically reduce to the well-known Routh-Hurwitz’s conditions, 
which we will extend to include not only asymptotic stability but also stability and 
therefore uniform stability. 


Hurwitz polynomials are intimately connceted with the problem of stability of 
mechanical systems and electrical networks and have numerous practical applications. 
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2. DEFINITIONS AND NOTATIONS 


Definition 2.1—The polynomial f (A) = A” + a, A") + ... + a,-,A + a, with 
complex coefficients is a Hurwitz polynomial if all its roots have negative real parts. 


Definition 2.2—If g (A) is any rational function, its paraconjugate is defined by 


g* (A) = g(— A) where A denotes the complex conjugate of A. 


When f (A) = A" + a, A"-2 4+ OW + a,-2 A*> + a,,A+ a, 
then 
F® (A) (= DEAR 1) ati A ee d,-, A* = a, , A+ a. 


Definition 2.3—A function h () is said to be positive if Re h (A) > O whenever 
Re A > 0. 


Theorem 5.! of Levinson and Redheffer’ reduces the study of Hurwitz polyno- 
mials to the study of positive functions, since it states that if f and f* do not vanish 
simultaneously, then f is a Hurwitz polynomial if and only if the function 
h(A) = joa is positive. This function Satisfies h* (A) = — 4 (A), since 
f** = fand hence — h* js positive if h is. Theorem 5.2 of Levinson and Redheffer’ 
States that any rational function h such that h and — A* are both Positive can be 
written in the form: 


b b, b 
tel eh py el 
— io 


A — Io 2 py _ iw, 


A()=atbr + 


where Rea = 0,b > 0, & > 0, and where the ®; are distinct real numbers. Accord- 
ing to the proof of this theorem, the ©; are the roots of 1/h (A) and since they are dis- 
tinct, it is an easy exercise to Show that the above expansion of h (A) is unique. The 


result of this theorem is a systematic decision procedure which will be heavily used in 
the next section. 


3. Asymptotic STABILITY WITH ComPLEx A 
We need the following : 


Lemma 3.1—If f(A) = An try ANSE ae ek a,_, A + a" is a Hurwitz polyno- 
mial, then f and f* cannot have a common root : 


PROOF : Write f (A) and f* (A) in the factored forms : 


and 


IA =(— I" ALAA & AD seen 


DIFFERENTIAL EQUATIONS 965 


Assume that A; = — a, is acommon root of f and f* for some i and j between 
1 and n. 
Then 4; + A, = 0 which leads to Re A; + Re A; = 0. 
This last relation means that A; and A j cannot both have negative real] parts, which con- 
tradicts the fact that fis a Hurwitz polynomial. 

Now, we consider the case where A is 2 x 2 complex matrix with characteristic 
polynomial 

f(A) = + aA +a, Let f* (A) = A? — 4, A +4 4, 

be the paraconjugate of /. 


Theorem 3.1—The system XY = A X where Ais a2 X 2 complex matrix with 
characteristic polynomial f (A) = A*® + a, A + a, is asymptotically stable if and only 
if: Rea, > Oand Rea,. Re (a, a,) — (Ima,)? > 0. 


Proor : X’ = A X is asymptotically stable if and only if f (A) is a Hurwitz poly- 
nomial (Jordan and Smith*, Theorem 9.3). Assume first that ¥’ = A X is asymptoti- 
cally stable, so by Lemma 3.1. f and f* cannot have a common root which leads to 
the fact that 

f(A) — f* A) 
h(a) = ST 
OSs Or are 


is a positive function (Levinson and Redheffer’?, Theorem 5.1). 
It is obvious that 4 (A) can be written in the form : 


(Re a,) A + i (Im a,) 


LI es A? + i (Im ay) A + Rea, ~ 


For any complex number « + 0, it is easily checked that Re « and Re (1 /«) have 


the same sign. Hence, 1/h is positive if h is positive. And since h* (A) = — hA(A), then 
(1/h)* = — 1h, hence — (1/h)* is positive. Therefore the function 
1  A®?+i(Ima)A +Rea 
h(a) (Re a,) A+i(Im — a.) . ’ 


can be written in the form indicated in Theorem 5.2 of Levinson and Redheffer’. 


We claim that Re a, ~ 0. For if Rea, = 9, then 


_ A? + i(Imai) A + Rea, 
a i (Im a,) 


E 


> 


(A 


where the degree of the numerator exceeds by 2 that of the denominator. 


~~ 


Therefore 1/h (A) can no more be written as in Theorem 5.2 of Levinson and Redheffer’ 
where the degree of the numerator exceeds that of the denominator by at most one. 


Hence Re a, # &. 
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By executing a long division, we bring 1/A (A) to the form: 


; 1 
| l | 
See eee ee I — Ima,) + 

h (A) (Re a:)® ig an Re a, 


Im — a 
+Re a.+ mer (Re a, Im a, — Im a,) 
(Re a) A + i (Jm a*) 


When we compare this form to that of Theorem 5.2 of Levinson and Redheffer? 
which is unique, we come to the conclusion that Re a, > Oand: 


Re a, (Re a, Rea, + Ima, Im a.) — (Im a)? > 0. 
The latter inequality can be written as : 
Re a,. Re (a, a,) — (Ima,)? > 0, 
We claim that Rea, Re (a, a.) — (Im a)? ¥ 0). 
for otherwise 1/h (A) will have the reduced form : 


l 


l 
ta Raat (Rea, Ima, —- Ima,) + 
1 


Re a, 





‘A 


which means that f — f* and Jf + f* have acommon root, which in turns implies that 


f and f* have a common root, contradicting the fact that fis a Hurwitz polynomial 
(Lemma 3.1). 


Conversely, if the inequalities of the Theorem are satisfied, 
I/h (A) as in the first part of the proof, we conclude 
Redheffer’,, Theorem 5.2), hence A is positive. It is enough to notice that f and 5 Ras 
cannot have a common root, for otherwise f — f* and f + J * will also have a com- 


mon root, which implies that the numerator and denominator of 1/h have a common 
factor, hence altering its form. 


then by writing 
that 1/A is positive, (Levinson and 


Hence fis a Hurwitz polynomial 


(Levinson and Redheffer?, Theorem 5.1) and the 
Proof is thus complete. t 


Now, consider the case where 4 is 


a3 x 3 complex matrix with characteristic 
polynomial 


F(A) = a8 + a, ae + a, A + as, 
and Let 

f* A) — A 4 3, rz a, A + us be the paraconjugate of f. 
Theorem 3.2—The system XY’ 


characteristic polynomial J (A) = A3 
and only it : 


= AX where 4 isa3 x 3 complex matrix with 
+ a, * + a, A+ ag, is asymptotically stable it 


DIFFERENTIAL EQUATIONS 967 
(1) Rea, >0 
(2) Rea. Re (a, a, — a3) — (Im a,)?> 0 
and 
(3) [Re a, Re (a, a,) — (Re as)? ]. [Re a,. Re (a, a, — a3) — (Im a,)*] 
— [Re a, Im (a, a3) + Re a; Im a,}*]) > o. 


PROOF : It is easy to see that 


h(a) = AL +i (Im a) 0° + (Re a.) + i (Ima,) 
(Re ai) A? + i (Im a,) A + Re az 


lige AG Ig asymptotically stable, then f and f* cannot havea common root 


by Lemma 3.1 implying that A is positive. Therefore h can be written as in Theorem 
5.2 of Levinson and Redheffer’. 


By a similar argument to that of the previous Theorem, we show that Re a; 0, 
After executing a long division : 


i I 
h(a) = (Re a,) (Re ai. Im a, — Im a,) + Reo 
x R 
(Re a,) A*+i (Im a,) A+Re az 


where 


] = 
R= (Re ai)? [Re a,. Re (a, 2, — a,) — (Ima,)* JA 


i : 
a (Re a,)? [Re a, Im (2, a,) + Re a3 Im a,], 


(1) is satisfied according to Theorem 5.2 of Levinson and Redheffer?’ The re- 
marks following this theorem as illustrated by Example 5.1 in Levinson and Redheffer’ 
allow us to extract (2) and (3) by repeating the argument of Theorem 3.1. The pro- 
cess is rather lengthy but straightforward and consists of another long division and 
one further application of Theorem 5.2 of Levinson and Redheffer. 


The converse is established along the same lines asin Theorem 3.1, we extend 
these results a little further. 


Theorem 3.3—The system X’ = A X where A isa 4 X 4 complex matrix with 
characteristic polynomial f(A) = A* + a, A® + a, A® + a3 + a, is asymptotically 
stable if and only if : 


(1) Rea, >0. 
(2) Rea,. Re (a, 4, — a;) — (Ima,)’ > 0 
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(3) [Rea,. Re (a, 4, — a,) — (Im a,)?]. [Re a). Re (a, 43 — a, 4%) 
— (Re as)*] — [Re ai. Im (4, a3 — a,) + Re as. Ima,]* > 0 


(4) a. b —c? > O where: 
a = [Re a,. Re (a, a, — a3) — (Im a,)?]. [Re a,. Re (a3 4,)— (Im a,)*] 
— [Re ay. Re (a, a,) — Im a, Im a,}?. 
=[Re a,. Re (ai 4, — a3) — (Im a;)’]. [Re a;. Re (a, a3 — ai 4) 
— (Re as)*] — [Re a. Im (a; a, — a,) + Re as. Im a,]? 
c = [Re 4). Re (ai a, — as) — (Im a,)*]. [Re a,. Im (a a,) — Re a3 Im a]. 
+ (Re a. Re (a, 44) — Imay Im ay]. [Re a, Im (4, a3 — as) 
+ Re a3. Im ay]. 


Proor : The procedure is entirely similar to the previous Theorems. A (A) can 
be written in the form : 


(Re a,) A® + i (Im a.) A® + (Re a3) A + i (Im a,) 


0) = aa a Uva) AP (Re ay) ae tT (meee AER 


In 1/h (A), the degree of the numerator exceeds that of the denominator by one. 
After performing a long division, the quotient produces (1) following Theorem 5,2 of 
Levinson and Redheffer’. In this division, the remainder over the divisor represents a 
second degree polynomial in A over a third degree polynomial in A, which when re- 
versed brings us to a position similar to that of Theorem 3.2 and the same argument 
applies. Again the procedure is lengthy and cumbersome, but straight forward. 


4. STABILITY wITH REAL A 


A well known Theorem by A. Hurwitz establishes necessary and sufficient condi- 
tions for the asymptotic stability of the system X’ = AY with real A. (Gopal’, Theorem 
8.5), Different approaches to the proof of this Theorem (sometimes referred to as the 
Routh-Hurwitz crierion) may be found in Bellman’, Cremer and Effertz* and Fuller‘. 


It is straightforward to verify how the results of the previous section, reduce to 
the Routh-Hurwitz’s conditions when A is real. In this section, we extend these condi- 


tions to include not only asymptotic stability but also stability and therefore uniform 
stability. 


Theorem 4.1—The system X’ = A X where Aisa 3 x 3 real matrix whose 


characteristic polynomial J (A) = 2 + a, a2 + a, A + a, has no zero repeated root is 
stable if and only if 


a, > 0, a, > 0, a; = 0, a1 a2 — az > 0 


where asymptotic stability occurs Only when all the inequalities are strict. 
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Proor : The case of asymptotic stability is settled by Theorem 8.5 of Gopal®. 
For stability which is not asymptotic, f(A) has one of the following forms : 
(a) f() =A (A — «) (A — 8) 
where « and 8 are real and negative or complex such that Re x < Oand Re 6< 0. 
(b) f (A) = (A — y) (A + 0) A — ia) 


where y is real and negative and @ a non-zero real number. (Jordan and Smith® 
Theorem 9.3) 


Case (a): \(A — a) (A — 8) = A® — (x + B)A2 + aBA 
= A® + a,dA2? +. a,rX + as. 
Hence, 
a + B = — a, « 8B = a, and as = 0. 


When « and 8 are real and negative, then a, > 0 and a, > 0, when they are 
complex with zero or negative real parts, then a; > 0 anda, > 0. 


Therefore if f (A) has form (a), then : 


a, > 0, a, > 0, a3 = O and these conditions imply that a, a, — a; > 0. The converse 
is easily established. 


Case (b) : (A — 7) (A + i) (A — 10) = A? — yA? + BA — Ye 
= M+a,A?+a,A + ay. 
Therefore a, = - yY, G2 = 6? and a3 = — y6* which leads to: 
a, > 0,a, > 0, a3 > Oanda, a, — a, = 9. 
The converse holds, and that ends the proof. 


It may be worthwhile for the sake of applications to state explicitely how the 
variations of the above conditions affect the nature of the eigenvalues of A. 


(1) a> 0, a,>0, a3 > 0, a, a, — a3 > 0. These are the conditions of nsymp- 
totic stability where the eigenvalues are real and negative or complex with 
negative real parts. 


(2) a, >0,a, > 0, a3 > 0, a, a, — as = 0, one of the eigenvalues is real and nega- 
tive, the two others are pure imaginary. 


(3) a, >0,a, > 0, a3 = 0, a; a, — a3 > 0. Zero isa non-repeated eigenvalue. The 


two others are either real and negative if a* — 4a, > Oor complex with negative 


real parts if aj — 4a, < 9. 
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(4) a, = 0, a2, > 0,a3 = 3, a, a2 — G3 = 0. Zero isan eigenvalue and the two 
others are pure imaginary. 

Theorem 4.2—The system X’ = A X¥ where A isa 4 x 4 real matrix whose cha- 
racteristic polynomial / (A) = A* + a, A? + a, dA» + a3 + a, has no zero repeated 
root is stable if and only if : 

a> 0, dyi> 0: as; > 0) as => 0 and a, a, a, — a; a, — a; > 0 
or, anein 
a, = 0, a2 > 0, a3 = 0, a, > O and as — 4a,>0 (2) 
where asymptotic stability occurs only when all the inequalities in (1) are strict. 


Proor : Again, Theorem 8.5 of Gopal! settles the necessary and sufficient con- 
ditions for asymptotic stability which can be simplified to appear as in the statement 
of our Theorem (Rao®, Remaks following Theorem 2.4 2). 

For stability which is not asymptotic, f (A) has one of the following forms : 


(a) f(A) = (A — a) (A—8) (A + iy) (A — iy) where « and 8 are real or complex such 
that Re « < Oand Re 8 < 0, and y a non-zero real number. 


(b) f(A) =A (A—a) A (A — 8) (A — y) where « is real and negative 6 and y are real 
and negative or complex such that Re B < 0 and Re y¥ =< 0. 


(c) fAJ=A—imwCAt inlA — i8) (A + iB) where « and 8 are non-zero real nu- 
mbers (Jordan and Smith*, Theorem 9.3) 
Case (a): (A — a) (A — B)(A + iy) (A — iy) 
= MM (@ + B) AM + (a8 + y°)A® — (@ + B) 2A + ay? 
= M + aA® + ae? + asd + ae 
Hence 


aj = — (a + B); a, = a8 + ¥?, a, = — (a +.B) y’, a, = aBy?, 


% and f are real and negative or complex with negative real parts if and only if 
“+B< Oandaf> 0, 


Simple manipulation leads to : 


a1 > 0, a> 0, a; > 0, a, > Oand a, a, as — a; dy Gag, 


It is easy to establish the converse. 


Case (b): d (A — a) (A — 8) (A — y) where « 


,B and y Satisfy the conditions al- 
ready stated. 
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It is obvious that aa = 0 for \ = 0 is an eigenvalue. 


Therefore f (A) can be written in the form : 
f(A) = A(A — a) (A — 8) — y) =A (A9 + a, + a2d + as). 
Theorem 4.1—Implies that a, > 0,a, > 0, a; > 0, and a, a, — a3 > 0 which 
when combined with a, = 0 leads to 


@, 4,43 — a; a, ~ a, > O2~ 


The converse is easily established 
Case (c) : f(A) = (A — ix) (A + ia) (A — 18) (A + 18) =A4 + (a? + B82) A? 
+ a? B® =)‘ + a, A? + a, A? + asd + a, 
therefore 
a, = a3 = 0, a, = a? + B*, a, = a? B2. 
«* and 6” are then the roots of the quadratic equation ¥* — a, ¥ + a, = 0 which has 


real positive roots if and only if a; — 4a, > 0, a,>0,a,>0. 


Therefore a, = 0, a. > 0, a, = 0, a, > Oanda} —4a,> 0. 
The converse also holds and that ends the proof. 


Again, we end this section by discussing how the variations of the above condi- 
tions (1) and (2) affect the nature of the eigenvalues of A. 


(1) a >0, a, > 0, a, >0, a4 > 0, aia,a3 —a, a,— ay > 0. 


These are the conditions of asymptotic stability where the eigenvalues are real and nega- 
tive or complex with negatixe real parts. 


2 2 


Here two eigenvalues are real and negative or complex with negative real parts and two 
are pure imaginary. 


(3) a, > 0, is > 070, > 0, @, = 0, a, a, a3 — a; ay = G, > 0. 


Zero is a non-repeared eigenvalue, the three others are such that one is real and nega- 
tive and two are either real and negative or complex with negative real parts. 


2 +R 
(4) a, > 0, a, > 0, as > 0, a, = 0, a1 4,0, — 2; & — a, = 0. 
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Zero is an eigenvalue, one real and negative, and two pure imaginary. 
(5) a =0,a, > 0, a, =0,a;> 0 anda, —4 a, = 0. 
All eigenvalues are pure imaginary. 
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The asymptotic behaviour as « ~ 0+ of solutions of the initial value problem 


dx/dt = u(t, x, y, €), x (a,¢) = A (€), 
« dy/dt = v(t, x, y, €), y (a, €) = B (@), 
Gta ft SD i=— 00 
where « > 0 is asmall parameter and x, y, u, v, A and B are doubly infinite 
dimensional vector functions, is established using the theory of first order 
differential inequalities in a Banach space. Under appropriate assumptions, 


asymptotic estimates for solutions of the above initial value problem are 
constructed in terms of solutions of the reduced problem 


3 x/dt = u(t, X, Y, 0), X (a) = A (0), 
v (t, X, Y, 0) = 0. 


The estimates so obtained contain boundary layer terms which explicitly 
describe the nature of the nonuniform behaviour of solutions as functions of 
tand «. 


1. INTRODUCTION 


The initial value problem (IVP) for the scalar differential equation (DE) 


ey” +f(tyy¥,9 =90,t 2 Jo: a (a, bj, 6 — = -(1.1) 
subject to the initial conditions (ICs) 
y (a, €) = A («), y’ (a, €) = B (€), eal Aca) 


wheree > 0 is a small parameter, was considered by Baxley’, Weinstein and Smith’, 
and others. These authors studied asymptotic behaviour of solutions and their deriva- 
tives of the above IVP as the small parameter € goes to zero, The results were extended 
to finite systems of first order ordinary DEs and higher order ordinary DEs by 
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Ramanujam and Srivastava®. Inthe present paper, using the theory of first order 
ordinary differential inequalities in a Banach space, we establish the asymptotic beha- 
viour of solutions of the following |VPase — 0: 


Q’ (X;¥) ey — vk, eye 
R’ (x,y): = (x(a4,9) = A(e) + O(¢) 
I, (2,<) = B(e) + O(€) (13) 


fo (x,y): =x -—u(4,x,y,)-= O1E J, 


where 
Sed Ce Sir rerio fo Crk re he ea 
Ws (- 5 May Vas Cig stds Veit ee eee 
KE (ty Kas 0) X50) ¥ =e Pete ee eee 
Am (5 Aa, Mos Aj. =)? B= Gee Boe i eee 


a prime’ denotes a differentiation with respect to t,¢ > Ois a small parameter, Z is 
the set of all integers, A (e) = A (0) + O (€) means A; (€) = A; (0) + O (e) with 0 as 
the Landau order symbol. The consideration of the above IVP (1.3) is motivated by 
problems arising when method of lines [4,5] is applied to singular perturbation problems 
for parabolic differential equations with a small parameter € multiplying the time deri- 
vative. Also, if one considers the DE (1.1) with the ICs (1.2) ina Banach space of all 


doubly infinite real sequences, one obtains IVPs of the form given (1.3) using proper 
change of variables. 


In section 2 we state and prove monotonicity theorems in a Banach space which 
Shall be used in the rest of the paper. Sections 3-4 deal with IVPs for second order 
linear systems defined in the Banach space and obtains explicit bounds for solutions 
and their derivatives. In turn, these estimates are used to discuss the asymptotic be- 


haviour of solutions and their derivatives. The results of section 3 are generalised to 
cover IVPs for nonlinear systems in section 5. 


2. PRELIMINARIES - 


In this section monotonicity theorems for first order ordinary DEs defined in a 
Banach space are developed. The later sections of this paper make use of these theorems 
in obtaining estimates and hence studying the asymptotic behaviour of solutions of the 
[VPs mentioned in section 1. The proofs of the present monotonicity theorems are dif- 


ferent from that of ones given in Walter‘. Also, they are generalization of the theorems 
presented for the ordinary DEs, 


Let B = 1© (Z, R) be the Banach 5 


pace of all doubly infinite real sequences 
with the norm of u € B given by 


lull = Sup Ju]. 
i€Z 
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If v,w € B, then v < wand v < w respectively mean 1 < wi,i € Zand v; < w; 
— 5,i € Z, for some positive real number 8. 


The continuity and differentiability of functions x (t) : J > B, can be defined as 


follows. The function x (t) is continuous at a point tf) € J (respectively differentiable 
with the derivative x’ (f,)) if 


Ix (to + h) — x (t))|| > Oash > 0 
(respectively ||(x (to + 4) — x (to))/n —x’ (to)|| > 0 as h > O). 


In the following C* (/), for a real interval 7 C R, stands for the set of all func- 
tions k times continuously differentiable in /. 


Definition 2.1— Consider a Banach valued function g : J) x D > B,D C B. Then 
the function g = g (t, z) is quasimonotone increasing in u if there exists a real positive 
constant M such that 


u<v==> g(t,u) — g(t,v) < M(v — u), u,v E€ B. ...(2 1) 
Theorem 2,2—Let u be a solution of the DE 

Pu: =eu — g(t,u,s) =0,t € Jo elec) 
that is, 

Pra «ou = 21 (t,.u, ¢) = 0,-1.6: 2 
subject to the IC 

Ru: =u(a,e)=A(),uEGU=C'(,) 1 CV), J: = [a, 5}. aAeet) 
Then for every v, w € U, the following implication is true : 


Pv <Pu<Py, thatis,P)v << Ppu< Piw 


...(2.4) 
Ryv<Ru<«Rw 
=> 
v(t,6<) < u(t,«) < w(t, 6), FE J, fash 2s5) 


provided that 
(i) g(t,z,6«): Jo xX DX (0, «,] > Bis quasimonotone increasing in z in the sense of 
Definition 2.1; 
(ii) there exists a real number 5 > O and a “test function” 
g(f,«): J x (0,<«,)—> B 
such that 
S; (1, €) = Sj (t, e), ixj 
si (t,€) > OonJ,i, 7 € Z, sEU 
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and 
ie (w+s,s)-Piw >4,5>0inJ, 06) 
Pv—Pv—as)>a,8>0,1E€ Z 


for every real positive number «, and for every but fixed «. Here ~,s = (..., % 5-1, a1 


Sep ee: 


(iii) there exists a positive constant L such that 
g(t,w + 8s,¢) — g(tiw + 18, €) < L s(t) (B — n) NRTA 
g(t, v — ns, €) — g(t, Vv — 6s, e) < L s(t) (8 — ») 20:3) 
for every real y, B such that 0 < 4 <€ B. 


Proor: In the following, a proof is given to establish the right inequality of 
(2.5). Similar analysis will yield the left inequality. If the right inequality of (2.5) is not 
true then there exist «, > 0,f* € J, andj € Z 


such that 
u, (t*, €o) — wy; (t*, €o) > 0. 
For this €, let 
P (t, €9) = inf [wy (t, €9) — ui (t, €9), 2 € it Gee 


which is a continuous function of J (Walters*, p. 98). Hence, there exists a point 4, € Jy 
such that 


Oe P (ty, €o)/ s(t, 98) a es [P il: €o)! 54580) ]>0, 
t 


where 


5 (t, €o) = se (t, Ey). 
We have 


u(t, €o) — [W (1, €9) + a 8 (t, €,) ] < 0. : she} 
Given » > 0 such that 


a ie aay] ts S(t, €o)(L + M + 38) «ual 2.10) 


and hence given ys (f», ¢,), there exists atleast one j/ € Z with the property that 


u; (to, €y) matty (fo, €o) —- as (to, €o) == nS (to, €o) noahieel 1) 
that is, 


Uj (to, &o) — [w, (fo, €o) + (a — n) 5 (to, €o)] > 0. ++o(2.12) 
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Hence the function 
uj (t, Eo) aa [w; (t, Ey) ss (a Man 1) § (t, €o)] 


which is negative at t = aand positive at t = f) attains a positive maximum say at 


Therefore at t = t, we have 


u; (ti, €o) — [W; (ti, €0) + (a — y) 8 (t,, €0)] > O Pel aes 
and 


ui’, (f1, €o) — [Wy (t, €o) + (a — y) 8’ (ts, €o)] > 0. (PH EY 


Now from (2.6), (2.14), (2.1), (2.7), (2.10) and (2.13) we have at t = 4 
0 > (Py w — Py (w + (@ — 4) 8))emr. + 5 (% — 2) 
eee ve Py Amat (% —y) 5), + 3 (2 — 9) 
= (vy; — w+ (qn — «) 5)’ (tr, €0) 
— (8; (t1, V(t, €0)s €0) 
— g,(t,w + (% — »)S, €&)) + S(a — n) 
> — (8; (4, V, &) — gj (t1,W¥ +S, Eo) 
+ 9; (t, wt *S, &) — By (tw + (a — 2), &)) + 3 (a — 9) 
> — M((w, + «5 — Vy) (th €0)) 
— Ls (th, &) +3( — 4) 


\/ 


2 — M(» S(t, €9))—L 1 S(t, €9) + & (a — 4) > 0: 
It is a contradiction and hence the proof of this theorem. 


The following theorem is stated for the system (13), 


Theorem 2.3— Let (x, y) be a solution of the IVP (1.3). Then for every x, y, y 
€ U the following implication is true 
[ P(x, y) <P (xy) <P’ Gy) 
Q’ (x, y) <Q’ (x y) <Q (x, y) 


Rwy) <R (wy) < RR Y) 


xixexyy Sy <y, 
provided that 
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(i) u(t, p,q, ©): Jo X B® x (0, €9] > B is quasimonotone increasing in p in the sense 
of Definition 2.1, and monotone increasing in q ; Y(t, p, q, €): 
J, x B® x (0, «] > Bis monotone increasing in p and quasimonotone increasing in q; 
(ii) there exists a positive number 5; and a “‘test function” s(t, «): J x (0,€] > B 
such that 
sj) =s(t,¢>01ij7E€ ZiFjsEUu 
and 
[Po + asy + a8) —P’ (x,y) > ud >0, 
| 
; at be cee alt ALi @,S,y — 4,8) > a, 4, > 0 


Q’ (x+a,s,y+ %s)+ Q’ (x, y)>% 8; >0 ..(2.15) 
LQ’ @, y) —O' @—“4, 8,7 — a,s)>a,5,>0 
for every positive real «, and for every but fixed e; 


(ili) there exists a positive constant L such that 

| u(t, x+Bs,y+6s,¢)—u(t, x+ ns, y+n5,«) 

: < Ls (t)(B — %), 

| u(t, x —ys,y — ns, €) — u(t, x — Bs,y — Bs, €) 

; < Ls (t) (B — 4). 

| (4x +Bs, y+Bs,c)—u(t, x +5, 

| y +s, «) <Ls(t)(6 — ») 

| v(t, x — 18s, y — ns,€)—u(t,x — Bs, 

l y —8s,e) <Ls(t)(8 — n) ...(2.16) 


for every n, 8 such that 0 < 1 <8. 


< 


Proor : Under the transformation Zo = Xt, Zsi41 = Yi, i € Z, the system (1.3) 
reduces to a single doubly infinite system. The hypotheses of the present theorem am- 
ount to the same assumptions of Theorem 2.2. Hence the proof of the present theorem. 


The following discussion for the linear system illustrates the above Theorem 2.3. 


Consider a doubly infinite system of linear second order ordinary DEs 


ex" +ax'+Bx=y,r1Eey, (Ze 
subject to the ICs 
x (4, «) = A (©), x’ (a, e) = B(e) ove( 2a) 


where 


a= (%;;) ay = 0, i fj, toe S Z 
ic ae (B,;), a x! = ( » 1 xy’, Koes x", or) 
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@x=(.., = 8B y 
yo x x y are Ls 
hoe ie lj pe. B, ; Xj, ) 


Y — ( - » Y—1s Yos YI, ae 
and a, B, y are assumed to be continuous functions in their arguments. 
Let y: = x’. Then the above system (2.17) — (2.18) may be written as 


Fe Gaya xe 0 

7 O,(x,y):=ey’ + ay+PBx=O0,1E J, 

| x (a, ¢) = A (e) 

| Ross: = | 

L y (4, «) = B (¢), ...(2.19) 
Theorem 2.4—The conditions (i — (iii) of Theorem 2.3 are satisfied for the linear 


IVP (2.19) provided that 


(i) B<O(ii) oy > d33>0,1€ Z 

see oo 

(iii) = By, >—LL>Oie Z. 
=—0O 


Proor : The assumptions (i) and (ii) of the present theorem yield the condition 
(i) of Theorem 2.3 whereas the inequalities (2. 16) follow from the condition (iii) of 
the present theorem. Finally using the assumptions (ii) and (iii) one can show that the 
function s, s; (t, «) = e*‘I* is a required test function, by a proper choise of k, for the 
system (2.19) satisfying the inequalities (2.15). 


3. ESTIMATES AND ASYMPTOTIC BEHAVIOUR OF SOLUTIONS OF LINBAR SysTemMS—I 
In this section we discuss the asymptotic behaviour of solutions and their deriva- 
tives of the IVP (2.17) — (2.18) by obtaining necessary estimates. Since the IVPs (2.17)- 
(2.18) and (2.19) are equivalent in the sense that every solution of one system is a solu- 
tion of the other and vice versa. We consider only the system (2.19) in the following 
study. 


In the rest of this section it is assumed that all the conditions (i)—(iii) of Theorem 
2.4 hold true for the IVP (2.19). Consequently the implication of the Theorem 2.3 is 
valid for the [VP (2.19). Further, in the following, M; =i = 1, 2,... are real posi- 
tive constants independent of «. 


Theorem 3.1—Consider the IVP (2.17) - (2.18) and assume that A (ce) = O = 
(ec). Then we have 


IIx (t, €)I] < llyll e"*-” hort) 
IIx’ (t, €) || < m |lyl| en"? Faloee) 
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for some positive constant m and for every solution x of (2.17)—(2.18). 
Proor: Consider the IVP (2.19) and let (x, y) beits solution. Define functions 
(x, y) and(x, y) as 
¥1 = |lyll ew", xi = — % 
Ji = m |lyl] e"'—9, yn = — J, 
where mm is a positive constant independent of « to be chosen suitably. 
We have 


Py (x, y): = ¥/ — J) = O = Py (x, y) 
oe ot ose) 2 
On (x,y): =e, + ani + pa, BUS 


i] , 
= (em? + aym + % Pri) Yi, emttes) 


IIvll > ye = Qu (x, y) 
by a proper choise of m. 
Also 


x (a,¢)>O=x(a0,y (a6) > 0 =y(ae), 
Hence, by Theorem 2.3, we have 


es y <-y. Bitton 


Similar procedure yields 


*SxX,yQy. 


...(3.4) 
Hence 
| x1 (t, €=) | < |hyl] eme-o ut -5) 
and 
| v(t, 2) | < m [ly] eme-a, ...(3.6) 


Taking supremum in (3.5)—(3.6), 


and using the fact that y == x’ we get the required 
estimates (3.1)—-(3.2), 
Theorem 3.2 Consider the IVP (2.17)—(2.18) and assume that y= O=A. 
Then we have ; 


Ix (t, €) || < «ll B (€) ||M, [emt-2) — e- FU a)1,] eai(3Z7) 


DIFFERENTIAL EQUATIONS IN A BANACH SPACES 981 
IIx’ (t, €) I< el] B (€)|] M, e™'-*) + |[B (€) |le~®,—- 91] ...(3.8) 
whers m > Oandm 5; — L > 0. 


Proor : Consider the IVP (2.19) and let (x, y) be its solution. Define functions | 
(x, y) and (x, y) as follows. 


¥ = € |B (e)I| [em—9) —e-F3¢-9)/,]/(m « + 83) 
Ji = ||B (€)|| [m € em) + 83 83" lg}/(m € + 8s) 


X= — £, 59 = — 31,1 € Z. 
We have 
Pu(x, y): =X; — ji = O = Pus (x, y) 
—— co 
BH YEO, or MPD By Fy 
=— 00 
> « ||B (e)|| e-9 [m & + mo, — L)/(me + 383) 
+ || Be) |] eye [— 82 + 82 + € L](me + 8) 
2 O = Qu (x, y), i € Z. 
Also 


x(a,«) =O < x(a, «) 


y (a,<) = B(e) < ||B(e) || = y(a, 8). 


Hence, by Theorem 2.3, we have 


<a SYA Y. ..(3.9) 
Similar analysis yields 
xaxX,y SY. ...(3.10) 
Hence we have 
| xi (t, €) | < el] Bell Mi [em — e780] (3:01) 
and 
(y(t, 2) | < € IB Ce) 1M, em! + IIB (€)I] e826" vas(3.12) 
The estimates (3.7) — (3.8) follow from (3.11)- (3. 12). 
Theorem 3.3. Let w bea solution of the IVP 


aw +Bw=y,w(a,e) =A (e) el oi os 
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with the properties that 


lw"l| < M, 
and 


lim w(t,¢o=ua(t), te J 
6-0 + 


where u is the solution of the IVP 
a (1, 0)u' + B (t,0)u = y,u (a) = A (0). 
If x is a solution of the IVP (2.17)—(2.18) then 
IIx (t, €) — u(t) I< ¢ M, em") 4 C (e) 
+ ¢ M, ||B (e) — w’ (a)|j [em*-) —e-8, 0-4), 
IIx’ (t, €-) — uw’ (t)l < ¢ Ms e™’-9) + D (e) + € M, ||B (6) 
— uw’ (a)l| em) + |B (6) — u' (a)I| e-83°"-9/,, t © D 
where C (e), D (e) approach zero ase > 0 + 


Consequently 
lim x(t,e)=u(t)reJ 
e>0+ 

and 
lim x’ (t,<-)=u'(t),t € Jy. 
e>0+ 

PROOF : 


Let v be the solution of the IVP 
evwt+av+BPr=y+ ew te J, 
V (a, €) = A (¢); v’ (a,¢) = B(e), 
From Theorem 3.1 we have 
IIx (t, <-) — v(t, €) ||I< « M, em(t-e) 
IIx’ (t, €) — v' (1, &)l| < € Mg em(t-a), 


Also from Theorem 3.2 it follows that 


WG) — wit, ll <€ Mell BQ) — w (a, 8)I| X [e“O-9) — 080-=0/] 
...(3.24) 


llv' (t, «) — w’ (t, €) I< € M,|| B (e) — w’ (a, €) |le™(-a) 
+ |B (€) — w’ (a, €) fle 854 -9/¢, 
Finally we have 
Ix (1, €) — u(t) Il < Ix (1, ) — w(t, ol 
+ [lw (t, €) — w(r)}I 


...(3.14) 


pAaak) 


(3.16) 


(a 


...(3.18) 


(3.19) 


...(3.20) 


..3.21) 


...G.22) 
...(3.23) 


Pa eh | 


.-. (3.26) 
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and 
IIx’ (r, €) — u’ (t)I| < IIx’ (t, €-) — w’ (1, ©) --.(3.27) 
+ |lw’ (¢, €-) — uw’ (t)]. 
The results (3.17)—(3.18) follow immediately from (3.22)—(3.27) and (3.15). 
4, 


ESTIMATES AND ASYMPTOTIC BEHAVIOUR OF SOLUTIONS OF LINEAR SYSTEMS—II 


In the previous section we obtained various estimates and in turn discussed the 
asymptotic behaviour of solutions of the linear system (2.19) under the assumption that 
the conditions (i)-(iii) of Theorem 2.4 hold good. The condition (i) of this theorem 
is a very important one and is imposed invariably for coupled systems of DEs. In fact, 
it is the quasimonotonicity condition for finite systems of DEs. Then if the condition 
(i) of Theorem 2.4 is not met by the system (2.19) one can still obtain the same limit- 


ing behaviour by adjoining a new system of DEs to (2.19) as follows. Following Miiller 
[8, 9] we adjoin the following system to (2.19). 


f — Sup[x’— p]=O,Inf [x’—p'] =O 
a, A 
| — Supley’ +ay+Bp]l=—y 
ne s 
: Inf [ey tay+Pyp]=y 
| = 
| — x(a,e)= —A(e), x(a,¢6) = A(e) 
L —y(@e=— BO), ¥@, =B 41) 
where 
H:={up'ix<eexuy <p’ < yh. 
Define 
+ feu if Bi; > O 
Ryo ; 
! O otherwise, 
De, oe Bi; = ie : 
Set 
A A A A 
Not = Ki, Xoler = — Xv Vat = Ibs Vets = — We 


Then the above system (4.1) may now be written as 
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( A ACA A A 
| 2A y:=x -y=O0 
| a A A A A oe) fi 
_ 
Qy0t41 (X, y) > = € Votqi + Oi Vein. — 2 B Xo, 
| : ja-00 4 
| 
OG! 3 A 
7 x B, Xo j41 Sag fF 
| ju—Cco 
| 
J A A Ag A A co - A 
Qii (KY) 2 = € Vat + Ot Yor — x0 BY Xaj41 


co ea AN 
+ 2 By, Xay = v0 
j~——0O0 


A 


Xot4 (a, €) a Ps A; (e), Xe (a, e) == A (€) 


| 
| 
| 
| 
| 
| 
| 
2 
7 : 
| Yates (@, €) = — Bi (6), Yu (@, 2) = Bi(@)i € Z. ..(4.2) 
It can be easily seen that the the above system (4.2) satisfy a condition similar to 


A A 
that of (i) of Theorem 2.4. Also if (x, y) is a solution of the system (2.19) then (x, y) 
defined below is a Solution of the new system (4.2): 


A A A A 
Xotey = — Xi, Xot = Xt, Voinn = Vi, Vo = Wnt E fA. 


The following theorem correspond to Theorem 3.3 in the unrestricted case 
discussed above. It is to be noted that the Theorem 3.3 was proved under the validity 
of the assumptions of Theorem 2.4. 


Theorem 4.1—Assume that all the hypotheses of Theorem 3.3 except the hypo- 
thesis. (i) of theorem 2.4 are satisfied for the IVP (2.17)—(2.18). Then also the conclu- 
sions (3.19)—(3.20) for the [VP hold true. 


Proor —Consider the IVP (2.19). If this system does not Satisfy the condition 
(i) of Theorem 2.4 then adjoin the system (4.2) to it. Following the procedure given in 
Theorem 3,3. We may arrive at the following estimates for the system (4.2) 


A A A EK 

| x1 — ui | <eMy e™'-*) + € M,l| B (e) —w’ (a)| 
[emt r) ae e-830—4)/,) + — (e) 

A A A A 

| Yor — Uy | < € My emt-a) + ¢ Mol B (€)— w’ (a)|| e™'-a) 
A 
+ || B(e) — w’ (a)]| e780 =, 4+ OE (e) 

A A A A 
Voter — Wary | € Mo e™-@) + € Moll B (6) — w’ (a)]) e('-o 


A A 
+ IIB (¢) —w (a) e839), + FICE Z 


where D (e), E (e) F(e) tend to zeroase > 0 + and 7 is the solution of the IVP 
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: . 
Vv ni A 2 A 

aja (t, O) usin, — Ss Br, (t, 0) wey + > By, (t, 0) uayi. = — VE 

| i=—co j=—0O 

< ed co 

| . zy s < A 

| %(t, 0) uy — S By (z, 9) ayat > BY, (t, 0) uy = ve 

| j=m=co j=m—0O 

| A A 

t u (a) = A (0), 7 € Z. (4.3) 


It can be verified that if (x, y) and u are respectively solutions of the system (2.9) and 


A A A 
(3.16) then (x, y) and u defined by 


A A Cis a 

Noi = Xiy Xoigs = — Mt Vat = Wty Voter = — Vi 
A A 

Ust = Uj, Usigt = — Uj, t = Z 


are respectively solutions of the IVPs (4.2) and (4.3). These observations complete the 
proof of this theorem. 


5. ESTIMATES AND ASYMPTOTIC BEHAVIOUR OF SOLUTIONS OF NONLINEAR SYSTEMS 


In Section 4 different estimates are obtained for solutions and their derivatives 
for the linear system (2.17)—(2.18). Guided by the experience with the linear system 
we now proceed to consider the nonlinear system (1.3) and derive estimates for solu- 
tions of the same. The estimates thus derived may contain boundary layer terms which 
explicitly describe the nature of the nonuniform behaviour of solutions of the nonlinear 


system (1.3) in ¢ and «. 


The reduced problem for the ! VP is given by 
- p’ — u(t, p,q 0) = O,t Eo 


v(t, p,q, 9) = O 
p (a) = A (0). (5.1) 


We shall now make the following assumptions. 
(ii The functions u, v of (1.3) satisfy the condition (i) of Theorem 2.3. 
(ii) There exists a function Y such that 
v(t, p, Y¥, O) =O Ae) 
and the resulting IVP 


FR oan, 
p (a) = A (9), 


(5.3) 
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has a unique solution X such that 


|¥' | <M ies 
fe a 
v(t, X—a, Y— 8, «)— v(t, X —a, Y,.d SKB re & 


where k > 0, for every nonnegative functions a, B. 


(iii) For some constant 1>0 and for every nonnegative functions a, a*, B, B*, a <a*, 
B < p* 
u(t, X + a*, Y + B*,«)— u(t, X +a, Y + B, e) 
— B) 


< 1(@* —a + B* 
0, Rae 2 6a ee ee 
< 1(a*—a-+ B* — B) ».(5.6) 
v(,X +8, Y + BY )—v(t,X +a, ¥ + B,¢) 
< 1 (a*—a + B* — 8B) 
v(t,X —a,Y Be pr SUES i ¥2— Gera) 
< 1(a* —a + B* — 8B). <at5.7) 


Theorem 5.1—Let (x, y) and (X, Y) be respectively solutions of the IVPs (1.3) 
and (5.1). Further assume that 


lu (t, X, Y, <) — u(t, X, Y, Ol Ce M, (5.8) 
and 
liv (t, X, Y, ell < «My. ...(5.9) 
Then under the above assumptions 
(i) — (ii) we have 
IIx (t, «) — X (t) || < « Ms em¢-% 
+ € M, ||B (€) ¥ (a)|| 
x [emta) _p-kia)/y) é; ...(5.10) 
lly (4, —) — ¥ (t) || Ge M, em- 
+ ¢ Mg|| B («) — Y (a)|| x ema) 
+ I B®) — ¥ (@)]| e*-9),, m > 0, tS: i15 


PROOF: Because of the conditions (i)—(iii) stat 
ed above, Th : 
cable to the IVP (1.3). Also we have ecorem 2.3 is appli 


IA()—- A(O) || Ge M. 


Consider the vector functions (x, y) and (x, y) defined by 
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x=X+ny=Yty 
x= X—y,y=Y—y 


where 
Ni bas hy Op X= REY) Xd) 
n (t, €) = € Mg em(t-a) 4 ¢ M, ||B («) — Y (a)|| X 
[em(*—a) — @~ka)/,) -=t(os03) 
y (4, <) = € Ms em'-*) + € M,|| B (c) — Y (a)|| ema) 
+ ||B (e¢) — ¥ (a)]| e~ke-aye, ...(5.14) 
We have 


eo Geey ys mit (8) 4 — & (1, X + n V4 7,0) 


— um (t, X, Y, «)] — m (t, X, Y, «) 
a My 4 — 1 (4-4), 

from (5.8), (5.3) and (5.6) 
> € em-a)[— M, +m Mz — 1 M3 — | Ms] 
+ ¢€ || B(e) — ¥ (a)l| en" [m M, — 1 M, —1M,] 
+ ||B (e) — ¥ (a)|| e-*-9 [k M, + le M, — 1] 

by (5.13)—(5.14), i € Z. 
Choose M,, Ms, M, and m such that 


Ma 1s 0): (GAS) 
k M; — Ms — M, — 1 M3 > 0 ...(5.16) 
—1M,+kM,>0 (5.17) 


m M; — M,— 1M;—1M;>0 
and 


m M, 7 ] M, 7 1M, = 0. TA ert? 


Therefore we have, 
P: (x, y) > O=P; (x,y). i € Z 
that is, 


P’ (x,y) > O = P(x,y). 
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Also 
QO; (x,y) =«¥, (te 7’ (46) 
— [wy(t,X +n, Y + y,<)— w(t, X, ¥ + y, ¢) 
+ v(t, X,Y + y, 6) — w(t, X, Y, )] 
— y (t, X, Y, «) 
> —<M,— «MM, bay sen ee 
by (5.4), (5.9), (5.5) and (5.7) 
> € emia) [__ M, — M, + em M, — 1 M, + k Ms] 
+ ¢€||B (e) — Y (a)|| em"-) [e mM, — 1M, —k M,] 
+ || B(e) — Y (a) || e*" 9, [-k + le M, + k] 
by (5.13) and (5.14) i € Z. 
Using (5.15) — (5.18) we have 


QO’, (x, y) > 0=Q, (x,y), i EZ. 


Hence 

Q(x, y > V(x, y). 
Again 

¥i (a, €-) = Xi; (a) + « Ms = Ai (0) + € Mg > As (€) = x1 (a, 6), i € Z. 
Similarly, 


31 (a, €) > yi (a, €), 1 € Z. 


That is we established the following inequalities 


P’ (x, y) < P’ (x, y) 
| 
Q'«&, y) & Q -(x, ») ' 
Spe | 
R’ (x, y) < R'(x, y) J (5.19) 
which, by Theorem 2.3, yield 
x(te)< x(te,y (te) < y(t,),1E J. .».(5.20) 


Similar steps yield 


X(N x(HO)y(t,— <y(t6,tEV/. 


The inequalities (5.10) — (5.11) follow from (5.20) — (5.21). Thus we have 
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lim x(t,«)= X(t), rE J 


e>0+ 
lim y(t,«) = Y(t), t EJo. 


e>0+ 


CONCLUSIONS 


The analysis of the previous section can be carried over to study 1VPs for higher 
order equations of the form 


ey = f(t, y, y™, ..., y°"), t © J 
subject to ICs 

y (4, <) = A’ (¢), ..., y"-”) (a, €-) = AM (), 
where « > 0 is a small parameter, 

y'") ; = d" yldt", Ai (ec) = A*(0) + 06), 1 E Z. 


The method adopted in discussing non-quasimonotone linear system can be used to dis- 
cuss nonquasimonotone nonlinear systems 
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For n 2 2, let P (n) denote the largest prime divisor of an integer. In this 
paper, we develop an elementary method for estimating the sum = 


251s x 
f (n) P (n) where f(n) isa multiplicative arithmetical function. 


1. INTRODUCTION 
Let P (n) denote the largest prime divisor of an integern > 2. Various sums 
involving P (n) have been studied by many, authors!-5.78 10-1215 Jp) this paper, we 
discuss sums of theform f (n) P (n) where f (n) is a multiplicative function. Our 
2<nSx 


methods are elementary and can also be used to estimate sums of the form F (n) 


maSX 
hods are 
it appears that elementary methods are 


& (n) where f (n) is multiplicative and & (n)is additive. Through analytic met 
available for estimating the latter Sums?’4?}3’14 | 


best suited for the sums discussed in this paper. 


2. PRELIMINARIES 
Let « (n) denote the Mébius function ¢ 


(n) be the Euler totient function, o (n) be 
the sum of all the positive divisors of n, w 


(n) be the number of distinct prime factors 
of n and fork > 1, dk (n) be the Piltz divisor function defined to be the number of 
ordered k-tuples (X15 Xa, s0.5-X¥) OF Positive integers such that X1, X2,X~, =n. Also let 
B (n) = eS p and B (n) = x xp where, as usual, P denotes a prime number 
and p*||n means that p*|ln and pei In. Further for any function g defined at primes, 
let the function G be defined by G(/) = Oand forn > | 


G(in)= = gp). 
“As ea) 


In the sequel, 5’ means : 
ASX 2<nSx 
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Lemma 2.1—As x — 09, we have 


Sy P(n)= > B (n) +0| << — | (2.2) 


nex nox 
and 
= Bin)= = Bn) + O (x log log x). -e-(2.3) 
nsx n<x 


Proor : We have 


he a> ae SothP LL, Bleep 
n<x pm<x,p(m)<p pmax pm<x 
p(m)<p,psv x p(m)< p,p)> Vx 
Bea On Dechy a Pay sp te aa D 
PMNxX,pNV x pm<x pm<x 
p<vVx 
ap i B (n) +o|> p= | 
n<x psvVx 


and (2.2) follows. On noting that 


= tai(ny ==" DB (n) 4 x ap 
nsx n<x p°m<x,~P2 
=) SOR) Ox S alp**] 
n<x Pp Sx,4+72 


we get (2.3). 


Theorem 2.1—Let f (n) be an arithmetical function such that f (n) < g (n) where 
g (x) is positive valued and increasing on [l, co]. Then 


= f(n)B() = D f (n) P(n) + O(g (x) x*?/log x). (24) 
n<x n<x 

x f(n) Bin) = = f (n)B (n) + O (g (x) x log log x) ...(2.5) 
n<x nex 

x f(n) B(n) = is f (n) P(n) + O (g (x) x* "flog x). ...(2.6) 
n<x n<ax 


Proor : Since 
= f(nyBb(n)=  f(n)P (n) + & f(n) (8 (mn) — P(n)) + O()) 
n=x nx n<x 
= Bf) Po) + O(e@) Er) Pe)+ OM) 
n<x nsx 


(2.4) follows from (2.2). Similarly (2.5) follows from (2.3) and (2.6) follows from (2.4) 
and (2.5). 
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Remark 2.1 : Theorem 2.1 extends and sharpens an earlier result due to 
(De Koninck and Ivié®, Lemma 3). 


Lemma 2.2—Let G (n) be as defined in (2.1). Then 
k k 
ZL d(n)G(n)= =% (—1) ( i ) [ % g (P) d& (n)). 
n=x f=1 pin<x 


Proor: We ius Sf Ivic and Pomerance’®, p, 4-5) that dy (n) is multiplicative 
a — . 
and & (p*) = [ r for prime p and positive integral «. Also we note the 


following identity : 
ak k+j—1 
(pif ] [ a |- Ofork >landn>1. ...(2.7) 
i,j70,i+j=n 


This follows from 


ial 6 faced 3 as Obey 3010 {> Gua ; | a} {> eae | =} 
‘4 = meelek + j— 
Pe age blll 
Now by (2.1) 
ze tk (n) G (n) wate d,. (n) ee 8 (P) ts Pee (P) dy (pn) 
" ae ee Gh aes &(o) pA) 
-[*] > & (p) de (n) aie : (P) de (n) | 
Hea » & (p) dx (pn) 
-[1] > eam + > w {i (Pn) 
pn<x p’n<x 


-{#]a00} 
-[ va > eam +z, ..s(2.8) 


pnsx 
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say, Now since| “$' ]—[ 1 p= - isi} we have 
So he Ss & (p) dk (n) + », g (p) |e (pn) 


PN<x,P }n p?n<x 


a | | [acorn —=- {5} SS g (p) ad (n) + ») g (p) 


p’n<x p?n<x 


| (p*n) — i : ] dy (p*n) + [ : | dk (pn) \ 


Continuing this, we arrive, by (2.8), at 


» mom => coef *] >: g (p) dx (n) 


nsx i=] pin<x 


+ SS g (P)} di (pt)—| ; ] di (p*n) + += 5 Jaton 


p*, In&x 


However, this last sum is zero since for any prime p and any integer « z 0 


d:. Cols e) -[ : | dk. (p*t*) + aa + (—1)* | i |e Cre) — 0. 
In fact, this is equivalent to 
k+l+a+k—1]_ (k tiem a TF] 
[ k+1+« ] re k+o Beer Aoi k 
: l+a+k—-1 ] ey, 


1+ 


which is the casen = k+1+« of the identity (2.7). This completes the proof of the 
Lemma. 


Lemma 2.3—Let G be defined as in (2.1). Then 


= p(n)G(n) = — = 8 (p)e™) ..(2.9) 
n<x p*n<x 

= ¢(n) G(n) = = (p— Ig (p)¢ () ...(2.10) 
n>x p*n<x 


and if F is any multiplicative function satisfying 


F (p**2) + pF(p*) = F (p) F(p*™) aL) 
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for all primes and positive integers «, then 


x F(n)G(n)= = g(p)F(p)F(n) — = pg(p) F(n). ...(2.12) 
n<x pnsx p?n<x 


Proor : We prove (2.9), the proofs of (2.10) and (2.12) being similar. 
Zow(n)G(n)= S p(n) & g(p)= =X g(P)pu(pn) 
n<x p\n pnsx 


nx 


— y a as. 
pn<x,p fn & (Pp) # (n) peo & (Pp) » (n) 


en wae ee 
PrSx,p|n 


=— % g(p)u(n)+ = g(p)p(pn) 
pn<x p’n<x 


=— = g(p)u(n)— , =  &(P)e(n) 
pn<x PnSx,p} 


X,pyn 


— = g(pr@™— 5 g(p)u(n) 
pn<x p?n<Sx 


+ 2% g(p) p(n). 


9 
te x 


Continuing this, we get (2.9). 


3. Main RESULTS 


In the following, N denotes an arbitrary but fixed positive integer. 
2 


Theorem 3.1—There exist constants ao = Aa » @, ..., GN-; Such that 
N 
9 Qy- 9 
> PQ) =x? > Torr + 0 ("log x)" (3.1) 
nex I=] 
N . 
aie 
DS Bin) = DS Goer + O (x*M(log x") ...(3.2) 
nx i=] 
and 
N 
Ss B (n) = x8 oes aye + O (x*[(log x)¥*, 3.3) 


nx f=) 
Further, we have 


a = | tog ot 


1 


...(3.4) 
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where [x] denotes the largest integer < x and 
Sha eS Oy Se ay se ae 


Proor : We prove (3.2) and note that Lemma 2.1 and (3.2) yield (3.1) and (3.3). 
By the prime number theorem and partial summation, we have 


t x 
> = Sree + lage |) aCe 
PX | 
Hence by partial summation 


= B()= FS p+ TOC = Dl 
nx pn&x MSx/2 p<x/n 





x/n 


Pailin sealer 
na log t n® (log x/n)N*} 


n<x/2 2 





*/2 


ne [1] eel 
Paty i log (x/1) ab o| (log x)V*1 | 


1 
x/2 


N-} 
ee [t] log ¢ | | log i} ]} O Fes 
~ log | EDS [ log x 7 log i Line (logx)¥ +1 
a | | 


1 





oo 


war Sorel] MY a+ o [4] 


1 
x? 
cna | (log x)" ] 


i a; — | , x? ] 
-# San +l tora 
i=1 
where q, is as given in (3.4). 
To prove (3.5) we have by partial summation 


x 


[x] (t] 
Si. = aed 2 | 3 dt. 
nx . 

So that on letting x ~ co, we get ¢ (2) = 2a. By a similar argument, we find 


Co , 
53 MOS I ahs es fat SL veto?) 
- 


n=) 
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ioe) 
l 3 
It is easy to show that ¢ (2) < 2 ay ae so that dg < a. Fori = 2, we have by 


n=l 


(3.7) 


acd . 
2a; = iaj-, + >a = inte 20l-4 


n=2 


since a; => 0 for all i. This proves (3.5) and thus the theorem. 


Remark 3.1: Theorem 3 | is due to (De Koninck and Ivic®, Theorem 1) who 


sharpened an earlier result of Alladi and Erdés?. De Koninck and Ivic proved (3.1) 
and deduced (3.2) and (3.3). Our method of proof has an advantage in that it yields 


the representation (3.4) (and the inequalities given in (3.5)). 


Theorem 3.2—There exist constants dy, a, ..., dy_, such that 


N 
, of eee dy-, x? 
D>, de Pm =x >, tas + toa | 
i=1 


nsx 
Lente ; 
DORM : > ieee ol ear | 
gets : 
a d, (n) B(n) = x? 2 ine +O l aor | 


Proor : First we prove (3.9). By Lemma 2.2, we have 


y sme@-S eve [F] SY vue 


nx é p’n<x 
k k 
[I] DS moe Sem [4] 
pnSx t=2 , 
Xx >> pa, (n) 
pin<x 


=14 Sen fA]s 


ves( eS) 


roe ee E 


...(3.10) 


..3.11) 
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say. Since = dy, (n) <;. x (log x)*-) (cf. Ivic and Pomerance’’, p. 263), we have for 
x: 


n 
BS.2 


2 = 2 pdx(n) = =p = dln) 
saps pi<x — n<x/pi 
‘ Sy p( ~- ) (108 alle < x (log x)*-! log log x. (3.12) 
psx 


For the estimation of &,, we note that by (3.6), for each positive integer r, there exist 
constants 5,, b,, ..., b, such that 


: : b, x 
>= » Toe Cee eek 














px 
Hence 
X= FT pdk(n)+ = pd, (n) 
pnSx,nSvV/x PnSx,n>/x 
2 
a > d,, (n) » en ol > ds (n) | 
NV x p<x/n n>4/x 
x e x7, —_(x|n)? } 
= » fay n > bi log n ] ee Ee x/n)t} 
nv i=] 
+ 0 (x12 (log x)*~) 
dx fa (n) _ _logn ie of x ] 
- Sater (log oda i> log x me (log x)r*} 
nS a4/ x 
(S4e(S [ee] 
-»> Tera 5; /|_ log x 
nNvV/x yu 
eel ltt Olas J 
se LLicg =a | ys (log x)'** 
— j k-1 
. by bij d;, (n) (log n)’ 0 (log x) |} 
= x (log uJ > n* x1 /2 
1<i<r i 
osjar 


x? 
"2 Aiea | 
= Cy ety itt Me 
z Saray 53 o| (log x)"*? 
i=) 
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Now (3.9) follows from (3.11), (3.12) and the above. 


On taking / (n) = d, (n), g(x) = x* and recalling that d (n) <s n* for each e>0) 
in Theorem 2.1, we obtain (3.8) and (3.10) from (3.9). 


Theorem 3.3—There exist constants e;, f, 21, 0 < i < N — 1 such that 














r 
TORO ee 2 aay + o| aor | ...(3.14) 
" 
ai Pe fee ra | aces | ...G.15) 
> p(n) Bin) = x? y Tide aa O | (sear | ...(3.16) 
N 
: é (n) P(n) = x3 2 aay + ol eaer | _.G.17) 
Ny 
eS é (n) 8 (n) = “St lor | ane 
r 
Ps EO) Bla aap 2 a + O| tog sy | 3.19) 
y 
PS: SPE ae ae + 0 Paes ] nen 
N 
DS eg eats Pe Tere + 9 beer ] 3.21) 
and 
r 
Ps (n) B(n) =x > iene apie: _ | (3.22) 


Proor : The proofs of (3.15), (3.18) and (3.21) are similar to that of (3.9). We 
use (2.9), (2.10) and (2.12) in turn, instead of Lemma 2.2. 


The remainin i 
of the theorem follow from Theorem 2.1. 2 aoe 
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4. OTHER RESULTS 


Our general summation formulae given in Lemmas 2.2 and 2.3 and our methods 
of Section 2 can be used to discuss sums of the form = ff (n) g (n) where f (n) is either 


415X 
dy or « (n) or ¢ (n) or o (n) and g (n) is an additive function. As a further example, we 
have 


Theorem 4.1—For any integers k > 2 and N > 1, there exist constants Aj, B; 
andC;,, 0<igk— 1and0 <j < N — 1 such that 


k-1 k-) 
x dx (n) » (n) = x loglog x & A; (log x)*-1- + x & B, (logx)*-$ 
n<x t= 0 7=0 


N-1 


+e Dora * Olea |. 


i=0 
ProorF : Jt is well known (cf. Ivié and Pomerance?®, p. 163) that for certain 
constants ae Oe a < e 


k-1 ; 
2 dg(n) =x 2% a? (log x)**-* OO [x?-*"], ack) 
nSx bead 


Now by Lemma 2.2 


Baio = 3 (9 | FY aon 


-| Ble D (5 13 eal S, (4.2) 


say. 


Fori > 2, we have by (4.1) 


k-1 
a 
Sao S (ese eh 


pxx'it n<x/p* p<xi/t 
x 1-1 fk \ 
ps | | 


+0f 
1(— kK Vey 
2 Se SAS cre own 


j= psx} it r 


M 
I 





r ae Ty i(1-11k 
x (log p) } + O[x a P ) 
(equation continued on p. 1000) 
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| | k-1-j 


Sot Sem osone ef ey 


j=0 r=0 psx'/* 
“be O.[x*o* log loge) 


hat k-1-j 





=>? Se [FTF ) coger {S log p) 
j>0 r=0 Dp 
zis of or x) }} + O [x1 log log x] 
od, Sp j (log x)*-2-> + O (x'4/* Jog log x] ...(4.3) 
s=0 
for certain constants D),0 <j < k — 1, 


For the estimation of 2; we use Dirichlet’s hyperbola method. 


x; + a dy (n) = D>» di, (n) + & dy. (n) — = d; (n) a x d' n) 
pnSx pnSx pnSx pn=x PSV Xn x 
n<v/x nV x 


=) ie eS yu (4.3) 


say. 
Now it is known by elementary methods that for every C >0 
| l 
nx 


where B is a constant. Hence 


Sa ae 5 dj, (n) = »: dx (n) 


PNSx, pV x PSVX nSx/p 

a x x 
“SAPS De ep solls} J 
PSV x — 

k-1-j : 
x st ory 

SDF Sel Jeon gon ty 

pRnv/x ro 


+O [ x40 ] 


(equation continued on p. 1001) 
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k-) k-1-j; 
= (k) | k-1-j r 
ep ey ee | Cor y eer 
7=0 j=0 


PSV x 


as O| xa | 


k-1 
= <x Das a”) (log x)*-1-/ {log log x —log2+ B+0O eal 

< (log x)¢ 

k-1 ae jo : 
man a'™ (— tap va] (lo k-1-j-y {98 xy 

54 x) rE + E, 
De aD: HA Eg 
] 
Bs Ol cere }} + O [xt-1/2h] .-. (4.6) 


Ks 
= x log logx 5 a (log x)" + x ; B® (log x)t-1-1 


3=0 j=0 


x 
eE o| (log pe*3*F =i 


For the estimations of and ie we write D, (x) = > d;,(n) and Ax (x) 


nwx 
k-1 
(k) k-1- 
= D, (x) — x 2 a, (log x)* =" 42 Then 
x? =| » ][ S ae ]=7Ww9 av 
pmv/x n</x 
[iv +0 [gals Tows 
= Li (Vx) De (Vx) + O [x/(log x] (4.7) 


=z 
0, by th ime number theorem. Here Li x -| aie Also 
for any c > 0, by the prim : fot 
2 


x/n 


oe = dx (n) S l =: dy (n) {Li | =] ihe ie 


; nxV/x pRx/n Nv x 


= S ds (n) (n) Lil x | a0 Feel (4.8) 


nsvVx 
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However, by partial summation 























/x 
: x , Ds (t) Dx 
d, Lil — ]=L D Pee ees 
> amul- ] = ties Dev) + BOF at 
nev x | 
Vx A 
=Li (4/x) Dy x | « (1) oes Ap dt 
x t log x 
1 
/x 5 
a® 
+x | aoa SST a, (log-t)*s*—7 i, 
1 j= » 
Now 
/ x / x 
x fe wa (t) _ _logt pa " | A, (t) gt 
Spee Be cg Pe ; Lets 
log x log x log x ‘ ! i ° 
1 





logt 4 log t |" 
Fi —— jN-} 
e aed + Olson | |e 


[ A, (t) (log t)! 


t? 


ee oa > “ons 





dt +- O (x (log x)¥). 


Since A; (x) = O [x'-1/*] Further on writing u for x/t, we get 


/ x 
k-y k 


= x 
x | —.—-) a (log 1)*!- dt = x >: at" | {log (x (u)}*7) 4 
9 


t log x/t S ae x/u log u 
x = du 
1 k-1- j : x 
==> a SU 7 7 dos ait | CEM 
7=0 r=0 Xx ms 
k-y kasd 


j=0 


: Sa | (log x)*~) (log 2) + Dy (— 1) ek wild ] 


k=] 
-l-r (I ad 2°") lo x) 
x (log x)* 1 i — 27) (log * } —— 9 5 F; (log x)*-1-3 


j=0 


for certain constants F,, 0 < j 
p> jJ&k-—1. Thus with a large C, from (4.2 
(4.8) and the above, we obtain the theorem. OS) ae 
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R . . 
oh sai : 4.1: The case kK = 2 of Theorem 4.1 was earlier attempted by De 
oninck and Mercier®. See also (De Koninck and Ivic?, Chapter 9). However, there 


. li . . 16 . 1) 4 . + 
a 


Ivic raised the problem of proving Theorem 4. | by an elementary method. 


We state the followin 

g results that can be proved by th 
ees y the methods of this paper. 
: of constants A, B, C, D and E£,, F,, ... » En-; may vary with the asymptotic 
ormula. For each positive integer N, we hits 


= 2-8 a 
met 2* (n) » (n) = Ax (log x) (log log x) + Bx log x + Cx log log x 


N-1 
eS ener pee 





D> emem=x gesy + Lge | 


nx i—2 


N-1 


> ¢ (n) » (n) = Ax? log log x + Bx? + x? > 
= (log x)’ 
i= 


nex 
OO ae ies 
+ Ol agg 7 | 


and 
N-1 
> 3 (n) w(n) = Ax? log log x + Bx’ + x* » ie 
(log x)§ 
ins i=1 


+ Olt a |: 


Some of these also appear in Ivic!’. Finally, we note that for estimating sums of the 


form > f(n)g(n) where f(n) is multiplicative and g (n) is additive, the method of 
nx 

this paper can be used when g (n) is an arbitrary additive function while it appears that 

the analytic method outlined in Ivic’* can be used only when g (”) is a prime-independ- 


ent additive function. 
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ON THE PROBLEM OF THREE GRAVITATING TRIAXIAL 
RIGID BODIES 
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In this paper we consider the problem of three triaxial rigid bodies formulated 
in Delaunay and Andoyer variables. Taking certain constant values for the 
energy corresponding to the approximation of three point masses, we can 
obtain the effect of interaction of bodies due to the non-sphericity on the 
Delaunay and Andoyer variables and hence in these cases the general motion 
of the bodies in the space can be determined. 


1. INTRODUCTION 


When we are dealing with triaxial rigid bodies instead of mass points, the 
problem becomes even more complicated. Aksenov'’ obtained solutions to the planar 
double averaged elliptical restricted three body problem. More recently Sidlichovsky®, 
considered the more general three body problem and obtained the solution of the 
equations of motion, assuming that the angular momentum of the close binary is much 
smaller than the angular momentum of the motion of the binary around a third body. 
Sidlichovsky*, obtained the Hamiltonian of the problem of three rigid bodies in 
Delaunay and Andoyer variables averaging over the fast variables, applying various 
restrictive assumptions obtained different approximations of the problem. Duboshin® 
showsthat the problem of translatory-rotatory motion of three rigid bodies admits some 
particular solutions (Lagrangian and Eulerian solutions) when each body possesses 
axial symmetry and equatorial symmetry. Cid and Elipe® studied the plane motion 
of three rigid bodies besides the case of three axisymmetric ellipsoids. In the present 
work we consider the gravitating of three triaxial rigid bodies M; (i = 0, 1, 2) such 
that the distance r of the bodies M, and M, is small compared to the distance r’ of the 
body M, from the centre of inertia of the bodies Mo, M,. Taking certain constant 
values for the energy corresponding to the approximation of three point masses, we 
can obtain the effect of interaction of bodies due to the non-sphericity on the Delaunay 
and Andoyer variables. 


2. SYSTEM OF THREE RIGID BODIES 
Let us describe the system of three triaxial rigid bodies Mo, Mi, M, with masses 
Mo, Mm, m, by Jacobi coordinates of their centres of inertia 7), 71, Tz. Let T be the 
center of inertia of the bodies Mo, /,. Let x, y, z be the coordinates of 7, in the 
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coordinate system Sp with the origin in 7, and axes parallel to the inertia system: 
Similarly, x’, y’, z’ are the coordinates of 7, in coordinate system S with the origin at 
T and the axis again parallel to the coordinate system. The rotational motion of each 
body is described by Andoyer variables’ L;, G; Hi, l, gi, hi (i = 0, 1, 2). 


Introducing 
en Vhs oi (mo + m,) me 
mo + m,’ My + m, + m, 


the kinetic energy T of the system is obtained as 
T= 4 m(# + y? + 2%) + bm’ (#2 + yp? + 2%) 
2 L 
17 Bin? Je. (COS IE te - i$ | 
3 > | Aw + BUH )@; =o Lyi) He Sore mee gh 


k=0 





where A™, B and C) are the principal moments of intertia of the kth body. The 
potential energy consists of three terms 


U = Ug + Ups + Ups +»0(2) 


where U;, is the potential energy of the gravitational interaction of the ith and kth 


bodies. Only a part U' of (2) will be taken into account in the unperturbed 
Hamiltonian 


; m,m - (Mm, + m 


where r is the distance 7, T,,r’ isthe distance T T,, and f is the gravitational 
constant. Introducing Delaunary variables Z, Gy) Hele, he EPG foe a 


that their relation with the Kepler elements a, a’, e, e’, J, I’ are given in Sidlichovsky‘, 


Let us assume that the mass mM, iS so distance from the orbiting system M,, My, 
that its attraction may be approximation by that of'a mass point. The Hamiltonian 
H* of the problem after averaging over the fast variables 1,1’, Io, Li, 29 and g, Is 


given in Sidlichovsky*. The equation of motion can be written 


in the canonical 
form: 


d . oH* 
Fy hi 8, 0, 'g’, Wo i, gi, ta) ee ee 
dt (8 8 is Bi, Au) 0(L, GH, L', G', HE, GW) (3) 


* oH* 
de (Ei OH, LGU Ty Gey = ee 
dt ( ‘ : 0 6 th &, h, Ke g’, i’, hh, 8, hi) 


H* = H*¥ + HY (i = 0,1) 
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where H* , H* are unperturbed and perturbed Hamiltonian which defined in 
Sidlichovsky’. 


3. APPROXIMATION INTEGRALS OF THE EQUATIONS OF MOTION 


Now, to integrate the equations of motion we use the method of small parameter 
8 which depend on the difference between the moments of inertia. In the zero approxi- 
mation we take 8 = 0 (A“) = BO) = C) the Hamiltonian can be written in the form 


2 p38 "2 '3 
pm p= m 
Beer ors. tte 


where R,» (Sidlichovsky‘), is the part corresponding tothe approximation of three 
point masses. In the zero approximation the elements of motion depend on the energy 
constant R which are given in Sidlichovsky* we shall study three cases in the first 
approximation : 


(a) R= R, = 3 [6S + 5r — (135rS)'"], S & 3/5 r 
G? = (5/3 rS)'!? = const. 

(b) R=R;,S > 3/Sr 

(c) R= Re =F GS — sr) S <r 
G? = S = const. 

where 

S = (G' — C)?,G’ > G 
r= L?, H+ H’ =C = const. 


In the first approximation (B # 0) the angles h, h', hi and their conjugates are 
the only quantities which may be possibly affected by the non-sphericity of the 
bodies. The perturbation of these elements can be obtained from the following 





equations : 
‘  0H;* . ,0H;* 
ae ——} dt 
w= | Craw Ge), 
by ty 
t aHi* ¢ (2Ai* 
(4 
Sh = tae ) ar, sH = — (s-) dt (4) 
- 9 
i* 
‘ 0H.* / 0H; 
pepe sm, -{ (Ja 
ee | ( ob; } ah bhi 
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where the subscript 0 refers to the values calculated in the zero approximation. 


We shall obtain the perturbation (4) for the cases (a), (b), (c) as explicit functions 
of the time t. The expressions obtained are given by : 


a 


3 pt m7? ICH eat A) — BOO EN rk 5 
n= Sh ee (feta) (2B) 
o 
i=0 m 6 G; Gi 


x {sin’ J, sin? 1° cos 2 [hig —- « (t — t)] — sin 2 I, sin 2/9 





cos [hi — « (t — t9)]} 245) 
sH’ = — 3 ps mm? SS (20 — 4()— BW) es 
‘< -)30 aa geen <_< es ) 
: (tig em) LG? Ss : 


Gi, 
x {sin*? J, sin? I, cos 2 [hig — a (t — to)] + sin 2/, sin 215 


x cos [hio — a (t — f)]} 


a7.(6) 








32 


=o eee § 
a 


ee gee BIKE, mT sin 21° 
ine (22H Ae = BO) ( a . Z a 
Ge (mp) + m) Li Ge 


‘sit 21. 0 
Fr Eas sin 2/; cos (Ajp — « (t — ty)) 
m Ly Gy 


u’4 m’7 sin Ay fe 


- 





vem’ sin? J, . , 
ous ay * sin® Tig cos 2 (hig—ax (1) 
(m, + m,) Ly G; m Ly G, 





BAS 
and “ 
P 3 ptm? 3h 49 
3 hoe > 4 22 40 gy (— ay ) 
3 9 
<0 om, L? G?7G?— FR 


G 
{[2sin 2/, (3 cos? lig — 1)] t + fl (sin 2/y sin? M4 
a 


Sin 2 (Aig — « (t—t,)) — 2cos 2/, sin 2/9 sin (hig — % (t—t))]} 
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e a (2C( — 4) — B® 
(m, + m,) Ls G4] Go 2s He ins ) 


0 


oF ip 


40 
— —— — i / 
( ot ] ) {[2sin 27) (3 cos? Tig — 1)]) ¢ + - 
10 


[sin 27) sin? Jip sin 2 (hig — « (t — fto)) 


+ 2cos 2/) sin 219 sin (hig — « (t — t,))]} 





22 OCH A pe) cD be 
Bhs a ( a ) { asin ae E m’ (3cos* J,—I) 
Al Gi, sd Hie Gig mi Ly G 
u'* m'7 (cos? I’ — 1) oe 
ee eck sin’ Tet sin? J) ain. -Ii9 


(m, +m,) L3 GP 
x sin 2 (hig — a(t — to)) + 2 (sin 27, — sin 2/,) cos 2; 
x sin (Aig — « (t + £))] \ 
where « be determined from zero approximation, for the cases (5), (c) equal to 


3 fm, ma*(1 — e;*) *F 
16 a aI G2 — Hy 


and for the case (a) equal to 


(sin 27) cos 27; + sin 2/J, cos 21,) 





3 im ma® (1 — e?) 2 (1 — 5 ey cos 280) 


] 13 5 9 
a i GA— HH? 





(sin 2/, cos 27; + sin 21, cos 21s). 


Thus, the perturbed inclinations of the orbital planes with respectt o the reference 
e according to (5), (6) and the perturbed orientation of the rotational angu- 
(i = 0, 1) in space can be determined from (7). The results show that 
from the perturbed Hamiltonian due to the non-sphericity of the 


plane chang 
lar momentum G, 
there is no effect 
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bodies on the eccentricities e, e’ and the semimajor axes a, a’ of orbits M;, M,. Also, 
the rotational angular momentum G; and its orientation in the body M; (i = 0, 1) 
remain unaltered. 
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In this paper, existence, uniqueness and continuous dependence of solutions 
to direct and inverse problems associated with a particular Cauchy problem 
for an ordinary differential equation are investigated. Such problems arise in 
the study of the water movement about a gas-storage reservoir of aquifer type 
when the corresponding influence functions degenerate. 


1. INTRODUCTION 


Consider a pair of real-valued finite functions p (t) and v (t) defined on the closed 
interval 7: = [0, 1] such that the intermediate function 


Q(t): = v (t)/p (1) in I, Qo: = vo/po SI) 
satisfies the initial value problem 

c 2a = p(t) — pp in/,Q(0) =), c > 0. mA 
Let 

P (0) = Po, p(t) > Oin 1 (3) 
and 

v(0) =, v(t) > Oin/ ...(4) 


where c, Pp and vy) are assumed to be fixed positive values throughout this paper. 
Moreover, denote by 


Dp: = {p © C(O, 1]: p(t) fullfil (3)}, D,: = {v € ClO, 1]: ¥() fulfil (4)} and 

b: : = {Q € C' [0, 1]: 2 (0) = Q)} the subsets of admissible elements considered 
WRT mes 3 ° 

below. Often it will be useful to prescribe upper and lower bounds on v: 


0 < Yuin S(t) S "max < 0 in /. ea) 


Therefore, we additionally introduce the set 
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D : = {v € D,: v(t) fulfil (5)}.. Throughout this paper, ||p|| : = ets | p(t) | and 
pi t 
|P||: = max ||P»|| denote the norms of an element p € C [0,1] and of a bounded 
Nyll=1 


vil= 
linear operator P : C [0, 1] ~ C [0, 1], respectively. 


If the water movement about gas-storage reservoirs of aquifer type is modelled by 
the influence function method*, the time-varying dependence between reservoir 
pressure p and volume of gas v may be expressed by the Volterra-Stieltjes convolution 
integral equation 


t 


| Soy dP) = peeve: ..-(6) 





0 


Here, the gas-occupied pore volume Q satisfies eqn. (1). The continuous monotonically 
increasing influence function F(t), with F(0) = 0 and F(1)=c> 0, possesses a 
continuous differentiable convex derivative dF (t) | dt < 0 whenever t >0 and expresses 
the reservoir pressure response to a unit rate of water influx. An analysis of direct and 
inverse problems associated with eqns. (1) and (6) is given by the author* for the case 
that dF (t)/dt € C'[0, 1]. This paper intends to give an overview of the extremal 
situation characterized by degenerating the influence function into a Heaviside form 


0 if ft=0 


This case is just described by eqns. (1) and (2). 


2. FORMULATION OF PROBLEMS 


There are four tasks according to the initial value problem (1)-(2) which will be 
of our interest in the sequel : 


(P1) — Find pE D, when c > 0, py, > 0, Ymax = Vmin > 0 and ve D, are 
given ! . 
(P2) — Find vy € D, whenc > 0, % > Oand p € D, are given ! 


(P3) — Find c > 0 when y, > 0,v,: = v(l) and p € D,are given ! 


(P4) — Find vy, > 0 when c > 0, Po > 0, p(t) > 
t € [t,, t.] and w(t): 
given ! 


O for a certain sub-interval 
=v(t)—v int € [yy fal 0 entire te ace tence 


The prediction problem (P1) is of direct n 
the identification problems (P3) and (P4 
Indeed, deviations of p (t) from the initi 
v(t), O < + <4 from the initial value v 


ature, whereas the control problem (P2) and 
) belong to the class of inverse problems?, 


al value py, are caused by deviations of values 
0° 
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We are going to study the existence, uniqueness and stability of solutions to the 
problems introduced above, where a particular a priori bound principle expressed by 
Lemma 2.1 is essentially helpful. 


Lemma 2.1—Let for given Vmax > Vmip > 0,p € Dy, andy € D, satisfy the 
eqns. (1) and (2). Then we have 








O< Oinis [> rin Qo < Q (t) < Siaay i “eas Qo in J - (8) 
0 0 
and 
in Vm . 
0 < Pain: = @™ po < p(t) < pmax: = “™ pp in I. 9) 
Vmax YVmin 


Proor: If p € D, and v € D, satisfy eqn. (2), then it follows © € Da. Now 
let the function Q (t) attain its absolute maximum over / at the point t pal ’ This 
implies dQ (t)/dt > 0 and thus p(t) > po. Consequently, Q(t) <Q (t) = v (t)/p (t) 
< Vinax/Po = (Vmax [¥o) QM =: Qmax in J and p(t)>¥ (t)/Q (t) 2 Vmin|Qmax = (Vmin / 
Vmax ) Po. If otherwise Q(t) <Q, in J, then Q(t) < Qnyax holds evidently. Inan 
analogous manner one obtains the remaining two bounds of formulae (8) and (9) 


A A i“ ; : 
from dQ (t)/dt < 0 when t> 0 denotes the minimum point of Q(t) over/. This 
proves the Lemma. 


Note that the bounds Qmin , Qmax 5 Pmin ANd Pmax are independent of the constant 


c > 0 in eqn. (2). Thus, the subsequent considerations can essentially be confined to 


the subset Dp: = {p € D,: p (t) fulfil (9)} of Dp. Let us complete this paragraph 
with the remark that due to Lemma 2.1 a steady state of the form v(t) =v, in J 
implies a steady state p(t) =p,in/. The converse is immediately obtained from 
eqn. (2). 


3. On THE WELL-POSEDNESS OF PROBLEM (P1) 
This study proves the problem (P1) to be uniquely solvable for all prescribed 


data and presents a stability theorem regarding perturbations of v € D,,. 


Theorem 3.1—For any given c > 0, po > 0, Vmax > Yimin = 0 andv € D,, the 
problem P(1) is uniquely solvable in p € Dy. Moreover, this solution belongs to the 


subset Dp. 
Proor : We consider the initial value problem 
23, Me ( ORS) ), pire git, 0:(0) =O, (10) 
on Fe NO (0) 


where the function v is extended to negative arguments ¢ < 0 by v(t) = %. sides 
the problem (10) is uniquely solvable in Q € C1 [0, to], fp > 0, then we have becaus 
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of Lemma 2.1 Qnin & Q (%). Now*set a: = 1 — ft, B:= Qyin /2 and 
2 


ee = ( 
Yo ems 





Vmax — + po). Obviously, 
a 1 (v(t) eel (“ us 
<+and 7 ("52 — po | Q, Po 





Quin 
(aw -) 








2, c 


4 V max 


rig 9 
aoe 
whenever | ft — t | & a, | Q, — Q(%)| < B,|Q, — Q(%) | <B and | Q(t) 
— Q(t) | < 8. By employing the Picard-Lindelof theorem® we thus obtain a uniquely 
determined continuously differentiable solution Q(t) in |t—t |< min (a, B/y). 
Since the quotient 8/y is independent of f, there is exactly one solutionQ € Dg 
satisfying the formulae (8) and (10). This implies the existence of a uniquely deter- 
mined positive solution p(t) = v(1)/O(t)in 7. Asa consequence of Lemma 2.1 this 


solution belongs to D,. This completes the proof. 


Because of Theorem 3.1 there exists a uniquely determined operator P: D, cic 
[0, 1] > D> C C [0,1] that transforms v € D, into the associated solution 
p:=Pve D, of problem (P1). 

Theorem 3,2—Let p'):= Pv) and p) := Py) denote the solutions of problem 


(P1) according to v") € D, and v® € D,, respectively. Then we obtain a Lipschitz 
condition of the form 








2 2 
P max P x 


Proor : We introduce the operator H: D, x D, C (C [0, 1})? + C [0, 1] by the 


‘ 


formula (H (p, v)) (1): = a _— mae | (P(t) — po) dtin J. There exist partial 
0 s 
0 


Fréchet derivatives 


1 


(0, H (p, ») p)(th= — () p(t) + p(t) dt 
p (t) 0 





and 
(8, H (p,) v) (t) = — Au: 
p(t) 


A A 
at any point (p,v) € D, X D, which are continuous ina neighbourhood of this 
point. It gets evident that the linear Volterra integral operator of the second kind 
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ep H (p,¥v) has a uniformly bounded inverse satisfying the inequality 
2 

A A Po x Bore 

I@> H(p, Wl <= (1+ —~) 


YVmin 





(12) 


A ‘ 


whenever p € D, andy € D,. Namely, all p(t) + | p(s) dt = q(t) in I is 














p? (t) é 
equivalent to 
do (t d 
AU = cid pit) = P(t)in/ ..(13) 
where 
e(t): = | p(x) dz, 2 = PO 50): = 0 
0 v (t) 
and 
p(t): Pw) 
+0 


We can write the solution of the initial value problem (13) by the explicit formula’ 


| 


9 (t) = exp (—a(t)) x J Bxoxate)) ma) de 


This implies 


— exp (-— att) <: Po x | exp (a(t)) 5 pecs) () q (t) dt PAO) q(t) 
y(t) ° v(t) v(t) 


p (t) 


and therefore the inequality (12). Moreover, 


\|d,, Hip. vil < es if p E D, aah ve De Hence, as a consequence of the implicit 
function theorem® we obtain for the Fréchet derivative P’ of P at the point y = By") 


+ (1-8) MED, 0< 6 <1, 


2 2 


A A A A Vv eae Pmax 
IP’ (v)Il <p H (Pv, ye HP» MIS se (! . ) 


Vmin 





Then the Lipschitz condition (11) results from the well-known remainder bound 
structure of Taylor Series®. This proves the theorem. 
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4. SOLUTION OF PRoBLEMs (P2) AND (P3) 


By integrating eqn. (2) one obtains 
‘ 
¢ (Q(t) — Qo) = § p(t) dz — py tin I. ...(14) 
0 


From this formula there are derived explicit representations of the solution to problem 
(P2), 


t 
v(t) = p(t) (2% —~ Pe ) + - | p(s) dz bins ..(15) 
Q 
and provided Q (1) 4 Q, of the solution to problem (P3), 


f p(t) dt — py 
Q(1)—Q, _ 


...(16) 


c= 


Hence, problem (P2) is uniquely solvable in y € D, if and only if the prescribed func- 
tion p € D, satisfies the inequality 


t 
f p()dt > pot —Q cini. cack DD 
0 


If otherwise (17) gets injured for some value t € I, then (P2) possesses no solution. 
Evidently, the prescribed values p(t) > 0 in the subinterval 0 St < Qo c/p, can 
never violate the inequality (17) since Po t — Qo c is negative in this case. That means, 
arbitrary positive continuous functions P(t)in0<t€&%&, may be generated by 
appropriate positive functions y(t)in0 <r< fo attaining the initial value v (0) = Vo 
whenever t > 0 is sufficiently small. 


Theorem 4,1—Let c > 0, ve > O and p € D, be given so that the inequality (17) 
is satisfied. Then, problem (P2) is uniquely solvable inv € D, and we have a local 
Lipschitz condition ; 


lL + Ip] 
GQ) — yl) Dae pr’ + Hel Q) — p(2) 
ae | (max | Hees - - lip P 
.-.(18) 
where, for fixed positive values c, Po and ¥g v") € D, and v®) € D, are solutions to 
problem (P2) according to prescribed data P € D, and p? € D,, respectively. 


Proor : The Lipschitz condition (18) follows from eqn. (15) using the triangle 
inequality. For all t € J, we have 


POO O1 <1 (a — 2 1) (p@ pm (Hp 


(equation continued on p. 1017) 
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t t 


ot = {p™ | p® (2) dt + p® (t) | p®) (t) dt 


0 0 


t 
S194 — 2 rpg — poo + 1 {omen | 
0 

+ 1p () ~ (=) | de + | pO) — pe | 


t 


[pm @ar< fia, — 1) 4 POTMOO 
c rs 
@ 


x | p(t) — p@(t) |. 


This implies the formula (19) and completes the proof. 


Note that the inverse problem 


(P2*) — Find vy € D, when c > Oand p € D,are given! is not uniquely solvable 
since ¥p cannot be determined by the given data. Indeed, if formula (15) yields 
positive values v(¢)in 7 fora certain value Q,, then v(t) also remain positive if Q, 
attains larger values of arbitrary magnitude. 


No we focus our attention on problem(P3). Apart from the steady state 
v (t) = vo, p (t) = Po in J the requirement © (1) ~ Qp is natural. Then the prescribed 
data determine a unique solution c > 0 whenever they are compatible with eqns. 
(1) and (2). 


Theorem 4.2—Let, for fixed vy > O and pp > 0, c?)? > Oand c®) > 0 denote 
the solutions of problem P(3) according to prescribed data p™ € Dp, v) (1) > 0 and 
p®) € Dp, v (1) > 0, respectively. Moreover, let 


AM: = vy (1)/p (1) — O, # 0 
and 

A): = v@) (1)/p@ (1) — Q, + 0. 
Then we have 


1 

vO) (1) v(?) (1) | (1) bs a) — p(2)]| 

oe yet SS p" (*) dt — po lp P 

my ~ pC 
G 

bce) — cI < AMT aa] SS ar Ol Gs 

(19) 





The triangle inequality is employed to establish the result of Theorem 4.2. Note 
that the estimation (19) reflects the various influence factors with respect to perturba- 
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tions in the data. For a stable recognition of the parameter c, first of all, the data 
must be chosen so that | A“ | and | A | are large enough. Evidently, A: = Q (1) 


1 

— Q, = 0 is equivalent to f p(t) dt = py, i.e., po is the mean value of p(t)in J. If 
0 

however v(t) > vp and v (t) 5& vp in [0, €] for all « > 0, then 2 (t) > Q, for all t > 0. 


In order to avoid small values of A, it suffices to force a sufficiently rapid growth of 
v (t) over the interval /. 


5. Some REMARKS ON PROBLEM (P4) 


A nontrivial task of practical significance is the stable recovering of the initial 
value ¥) when noisy data of p(t) and only increment values w(t) = v(t) — v, are 
available. We make some remarks on a special version of this class of problems 
expressed by (P4). 





Let 
; t 
(1) eg | xe ( = rl |) 
aM (ji= + El (pm) po) de, PO (i = AO 
0 
and 
+ l l 
N@) (2): = BT — (2) (x) — 1): aw 
rm BO + 2] OMG) — pe) de, pm @): OE 


0 
hold for all t € J and a certain value x > 0. Then we have 


v(t.) + « v (t,) (* (t,) + « v (t) ) 


PIE) — PMH) pM) ~ pO (a) 


‘ 
9 


I 
= 1 | (p(t) — pO) dt if0 S14 Sh <1. 


‘ 
1 


Hence, 


t . 
2 


(sar Pa. | ; 
“\p® (t,) — a | pate (p (1) — p® (t)) dt 


Sy 


> : 
+ +O) (Saray — pray ) 


] ] 
+0 )) | oe aa 
(VY + w (t.)) (5m (15) p®) (t,) ) > 
A Teen tare. ...(20) 
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bl — t; vie Ilw 7 
cg ASEH A (+ Bp team 
Vo l ] 
aa (t.)  — p@(1,) | 
Ilplle : = max | p(t) | t € [t,, te] av=(21) 


Theorem 5.1—Let exist a solution ») > 0 to problem (P4). Then this solution 
is unique whenever 


w(t) — w(t) se Ph) (P(t) — Po) 
Cc 


ary int. & [ta, tl. aed 24) 


PRrooF : Provided that p (t;) + p(t.) holds, the uniqueness of the solution to 
problem (P4) immediately follows from formula (20) by setting p(t) = p“ (1) = p® (t) 
int € [f,, t.]. Evidently, we obtain « = © and therefore vo > 0 as a uiquely deter- 
minded value from the prescribed data. On the other hand, we can argue in the same 
manner whenever there is at least one value ¢ satisfying t, < t < (t2) and pt; # p(t). 
However, this requirement is only injured if p (t,) = p (t) int € [t, t2]. This identity 
implies 


a() = 0(4) +2— — (p (4) — py) 


i. e., 





y= v (1) + P(t) (p (ti) — Po) (t — t,) 


v 
(1 : 


and 


ONE 2 ple neta EP et int € [f,, ty]. 


This provides the requirement (22) and completes the proof. 


Finally, note that formula (20) also yields the continuous depenednce of solutions 
to problem (P4) upon the prescribed data whenever these data are compatible with 
eqns. (1) - (2) and satisfy the assumptions of Theorem 5.1. 
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P-WAVE SCATTERING AT A COASTAL REGION IN A SHALLOW OCEAN 


P. S. DESHWAL AND NARINDER MOHAN 
Department of Mathematics, Maharshi Dayanand University, Rohtak 124001 
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The present problem can be idealised to study diffraction of compressional 
waves at a coastal region in a shallow ocean. The method of solution is the 
Wiener-Hopf technique. An exact solution is obtained in terms of Fourier 
integrals. The evaluation of the integrals along an appropriate contour in a 
complex plane results in reflected Rayleigh and scattered waves. Near the 
scatterer, the scattered field has the behaviour of a cylindrical wave and at 
distant points, it isa plane wave. The scattered field excites the particle 
motion at right angles to the free surface of the ocean and the horizontal 
component of the displacement is absent in the surface. Numerical computa- 
tions for the scattered field in terms of the distance from the scatterer and for 
the Rayleigh waves versus depth below the free surface reveal a sharp fall for 
the former and a slow decay in case of the latter. 


1. INTRODUCTION 


The problem of diffraction of elastic waves by a rigid plane barrier in a liquid 
medium has been investigated by various authors. The effect of a vertical plane 
barrier, fixed in an infinitely deep sea, on normally incident surface waves was first 
considered by Ursell' dealing with a two-dimensional case. Faulkner* has extended 
Ursell’s study to a three-dimensional case to consider the effect of a fixed vertical plane 
barrier of finite depth on obliquely incident surface waves. He has used the Wiener- 
Hopf technique® and has obtained approximate expressions for the transmission and 
reflection co-efficients. The problem of diffraction of compressional waves which strike 
obliquely at a vertical barrier of finite depth fixed in a liquid halfspace or in a liquid 
layer over a solid half-space has been discussed by Deshwal‘-*. Mann and Deshwa]? 
have studied the problem of Rayleigh wave scattering at the edge of a finite plane 
barrier in the surface of a shallow ocean over a solid halfspace using the technique of 


Wiener and Hopf. The scattered waves have the behaviour of reflected plane waves 
near the scatterer and of cylindrical waves at distant points, 


The present paper deals with the problem of diffraction of compressional waves 
at a coastal region in a shallow ocean. At a coastal region, there is shallow ocean 
over the solid mantle of the earth on one hand and there are hard rocks on the other 


hand. The shallow ocean is assumed to be a liquid layer and the hard rock to bea 


rigid plane boundary. The regionx <0,z>0 isa rigid medium in which there is 


P-WAVE SCATTERING 1021 


no movement of waves. In the region x < 0,z > 0, there is a shallow ocean in the 


form of a liquid layer of depth A occupying the region x > 0,z €h lying over a 
solid half space x > 0,z > h, 


2. STATEMENT OF THE PROBLEM 
The two media lie on the side of positive z-axis (pointing vertically downward). 
The liquid layer and the solid halfspace occupy the regions x > 0,0 < z < hand 
x > 0,z > A respectively as in the Fig. 1. A two-dimensional compressional wave is 


(0,0) 


LIQUID 





Z 


Fic. 1 Oceanic layer over solid mantle of the earth. 


incident at the rigid plane boundary and the corner scatters the waves. The wave 
equations in the liquid and the solid are 





Se EE pm 0, KH ky + ik, mi 
ox OZ. 

0° ad 49 / Ps ce 

Feet pe TR HO km hy + ik, (2) 


and a similar equation in). k andk’ are complex with their imaginary parts being 
small and positive. ¢, ¢: and 4, are the displacement potentials in the liquid and the 


solid media respectively. 


The incident potential satisfying the stress-free conditions is 


¢; = D sin (kz sin 0) exp (— ikx cos 8) wet) 
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such that the total potential in the liquid layer is 

be = bi + o (x, 2). (4) 
The displacement components (u,, 7) in the solid medium in terms of potentials are 
given by 


0¢, a Ov, od, 


a ort paki aif Pe i), 
“= Bx az? 2 : on; 


Since there is no transverse wave in the liquid, the v-potential is absent and the 
displacement components (u, w) in the liquid layer arein terms of the derivatives of 
¢ only. 


3. BOUNDARY CONDITIONS 


The conditions on the boundaries are : 


(i) ¢:(x,z) = 0, onz =O, forall x ---(6) 
ee Od; es re 7) 
(ii) aes 0; on x. =.0,:for0 << z <= h fa8 

(iii) bid — jaa! = 0,onx = 0, forz>h ...(8) 

Ox Zz 

: 0¢, Obr EN ~ Pol 
(iv) a + 5, = 0, onx = 0, forz Sh -.(9) 
(v) W = Wy, Pzz = (Paz), (px), = 0, onz =h .».(10) 


for x > 0, where the subscript | denotes the entities belonging to the solid and pzz, 
Pzx are the normal and shear stresses. 


We define the Fourier transform 


#, (Ps 2) = F 6 (x, 2) e* dx, p = a + 8 Regie 


and assume that, for given z, 
| ¢ | — D, exp (—k,| x |) as |x| > co het #3) 


where D, is a constant, then 
= Le ¢) 
| $+ (P,2z)| < J | & (x, z) | e-®* dx, welt) 


On account of the assumption, $4 (p,z) is bounded at infinity only when 8 > —ky. 


Therefore ¢, (P, z) is regular in the region 3 >—k, of the complex p-plane. die (p, z) 
and 1, (p, z) obey the same behaviour. 
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4. SOLUTION OF THE PROBLEM 
Fourier transformation of (1) leads to 


ft is 
ie (2) — VS (7,2) = (H) — pe 14) 
where 
Y= +(p? - kp, (15) 


The sign before the radical is such that Re (Y) > Oforalip. The subscript 0 denotes 
the value atx = 0. Let us change p to — pin (14) and add the new equation to it 
to find 


in [b4 (Pp, z) dy = P; z)] — y? [b+ (p, z) 4. ey (- P, z)] 


= D (2ik cos @) sin (kz sin 6). ...(16) 


The right hand side is the value of (0¢/@x)) from (7) and the complete solution of this 
differential equation is 


$,.(P, 2) + b+ (— Py 2) = Cy (p) e-¥? + 2 (p) e?? 
(2ik D cos 6) sin (kz sin 6) (17) 
2 p? — k* cos? @ 
Let us utilize the condition (6) to find c, (p) in terms of c, (p). The condition results 
in 


$, (p, 0) + $4(— p, 0) = 0. ...(18) 


Using (18) in (17) to have 
as — (2ik cos @) D sin (kz sin 6) 
$y (Pp, z) + $4 (= P, z) as 2c, (p) sin hYz — aap — k= cos? 6) . 
.-(19) 
Putting z = / in (19) and in its derivative with respect to z and eliminating ¢ (p) 
between the resulting equations, we find 


= = | =, =m) 
$y (p) + $4(—P) = > tanh Yh [6 (P) + $4 (—P)I 
(2ik cos 9) D sin (kh sin 8) 
Pa (p? — k? cos* @) 


(2ik? sin 6 cos @) D cos (kh sin 6) tanh Yh 
(p? - k* cos* @) ¥ 
(20) 
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We have used the notations ¢, (p), 4, (p) for ¢, (p, h) and ¢. (p, h). Similar 
notations will be used for the potentials in the solid medium. 


Let us take the Fourier transform of the conditions (10) to get 


= ik sin 0) D cos (kh si =, = 
$4, (p) + E Sn OD cos Chain ) _ 5. (P) — vs), 
p # kcos@ eee AD, 
+ err oat 

— Me [ $y (p) + PROB SRD) ae FO) 

+ 2 (p? ~ kb) $14 (p) — 2u, ip Ui, (p) ...(22) 

and 

— 2ip $,.(p) + (k’? — 2p") diy (p) = 0, onz = h, x > 0, - ...(23) 


A,, #1 are Lame’s constants of the solid and A is that of the liquid, k” being the wave 
number for transverse waves in the solid. Taking the Fourier transformation of (2) 
and using the boundary conditions (8) and (9), we find 
ee ~ = o 
dzi #1+ (P, 2) — $14 (—P, 2)] — YI [ds (P, 2) — Fra (— p, z)] = 0, 
foal ih St te) ...(24) 
whose solution is 


b14 (Ps 2) — $14 (— pz) = A (p) e+ Bp) er’, wiles) 


Since $14 (P, 2) — $14 (— p,z) is bounded when z tends to infinity, it follows that 
B(p) = 0 and 


bis (Ps 2) — brs (— Dy 2) = A(p) O. sat26} 
Similarly 
tie (Ps 2) — tha (—P, 2) = D(p)e 8 = + (pt — Re QD 


The criteria for signs of Y, and 5, are the same as for Y. Putting z = h in (26) and 
in its derivative with respect to z and eliminating A (Pp), we have 


mS = wih Te PY ne 
$i+ (P) — diy (— p) = Paw [os4 (pP)—-¢., (=p). ...(28) 


5. SOLUTION oF THE WIiENER-Hopr Typg EQUATION 


The functional equation (28) is solved by involvi 
ing th - 
We write (28) as y & the Wiener-Hopf technique. 
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Fic. 2 Contour of integration in lower half of the complex plane. 


¢,, (p) $4 (k’) ay ij 
Pet Sek See 
bis (P) T (yp kyP = aR Fees: aT yr aE | (=2k' (p ae k')'/?) 


4. (—P) 4 (k’) 
= $+ (- p) + ore Ge — kyr - ser | 


$i4 (k’) 


Qk (p= Ry Me 


The left hand member of (29) is analytic in the region 8 > —k, and the right hand 
member in 8 < k,. By analytic continuation, they represent an entire function. Each 
member tends to zero in its respective region of analyticity as |p| -—- °c. By an 
extension of Liouville theorem (Sec. 2.61, Titchmarsh®) the entire function is identi- 
cally zero, Thus we have 


$14 (p) = — 7 Fie (P) + $/, (k) 8 (p) ...(30) 


where 


I 


I 
g(?) = pe eye + (2k) te ash) 


Similarly 


Tie (=F Ga P+ Fin RDA) BD 
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where A (p) is obtained from g (p) by replacing k’ by k’. 

Using (30) and (32) in the conditions (21-23), eliminating $i, (p) and vy (p) 
between the resulting equations, $4, (p) + bs. (— p) is obtained in terms 4, (p) 
+ $, (—p). The result is used in (20) to obtain 


re =, ¢_p, — _ (2ik* sin @ cos 6) D cos (kh sin 6) 
$4. (p) “Hr ?, ( Pp) Pa p? a k2 Cos” 6 
2¢(p) (—A.ki + mw ¥?) $4, (k’) 


TO ...(33) 


where 














= tanh YA (2p* nl Sl Kost Ay kk’? + 2p) Yj) 4p* Pi 3, 
bh Ng ney aes ui MEY, Re OOF gE G4) 


From (19) and (33), it follows that 


is « sinh Y2 728 (P) (A KP +2 YF) bi, (KD 
$4 (p, z) + $+ (—P, 2) = F cosh Yi ay ae ea ey | 
(2ik D cos @) sin (kz sin @) 
POSS sp aarae Goe ga (22) 


The potential ¢ (x, z) is given by the inverse Fourier transform, i.e., 


oo +18 
l - 
¢ (x, 2) = =~ $1 (p, z) e-”* dp, x > 0 ...(36) 


—co +48 


where --k, < @ < ky. 


6. REFLECTED WAVES. 
The factor exp (—ipx) = exp (— iax) exp (Bx) in (36) makes the integral vanish 
8 - —oo in the lower part of the complex plane if x > 0. For waves in region 
x - 0, the contour is in the lower half of the complex plane where bi (—p, z) is 
analytic and hence 


co +iB 


1 = a 
¢(x,z) = Qn | [p+ (Pp, z) + $+ (—P, z)] ete dp, x > 0. -».(37) 
—0co +18 


The contour of integration excludes the poles p = + k cos @and includes the poles 
P = + Pn (n = 1, 2, 3,...) where -E Py, are the roots of the period equation f(p) = 0. 
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Schermann’® has shown that the period equation has a finite number of real roots 
P = + P,, n denoting the nth normal mode of the generalised Rayleigh waves, The 
conditions Re (Y) = 0, Re (¥:) = O and Re (5,) = O have been discussed by Ewing 
etal.!°to show that the branch cuts are hyperbolic. Im (Y), Im(¥,) and Im (5,) 
change signs along the two sides of the branch cuts as indicated in Fig. 2. The line 


integral (37) is taken along the real axis (B = 0) of the complex plane. Applying 
Cauchy residue theorem to the closed contour results in 
co B D 
j+j+J + JS + J +f = 2wi (sum of residues) ..(38) 
a eee Pec kee © 
where 
co —Pn-€} : —k cos §—e, _ _ kos 0-€, 
J = limit : + | + limit J + f+ limit J 
-° €,-0 -* Day €1,€2>0 —pater > €.70 —kcos 6+e, 
Pa— : ‘ ‘ 0o 
+ { + limit j a) mit f lo?) 
23 €},€2+0 kK cos 6+ €, zy €:>0 pate, 
€, and €, being the radii of circular indentations 2, Z,, 2, and 2, around p = — p,, 
P = Pn, P = —k cos 6 and p = k cos @ respectively. The contribution due to inde- 
ntations at p = + kcos @ are 
—+4 D sin (kz sin @) exp (—ik x cos @) ...(40) 
and 
+ D sin (kz sin 9) exp (ikx cos @). ...(41) 
And the contribution due to indentations at p = + p, are 
—ip x 
i(—A, k* + 2n, Y2,) 8 (pn) $i, (k’) sin Y, ze 
- yy euereremeet Ak" 791.) con Yehie ja 
n 
Yi, = [(p, — k)}'* .»(42) 
and 
res: : ee 
i(—A, k* + 2p, Y;\)g (—P,) $54 (k’) sin Y, z e 
3S mmr ka fi(-p,).coxYehio  . 
Y, = [(k? — P;. yr .-.(43) 


In (40), the wave is half the incident wave whereas (41) gives the wave reflected Si 
the rigid plane boundary and then from the free surface of the ocean. In (42) an 
(43), we have the generalized Rayleigh waves multiply reflected from the boundaries 
of the ocean. 
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7. SCATTERED WAVES 


Let us evaluate the integral along the branch cut Py. On Py, Re (Y) = 0 and 
Im (¥) changes signs along two sides of the branch cut. It is found that the integrals 
along two sides cancel each other. 


We now consider 


Ione) = 5 | B.D + Fe (— pe dp. (44) 


, 
Pr 


On px’, we use p = — k’ — iu such that Y¥, = +i (2k, u)'® = + i7,, K = 0. 
Integrating (44) along two sides of p,., it follows that — 


(x= 55 | UG (12) + Fe (0, My og — 18 (D2) 


+ $4 (—p, 2)y — yy ] ems * emu du. (45) 
1 1 
This is evaluated by using the result (Ewing et al.'°, p, 52) 
ry G (0) 1 (3/2) . G’(0) F(5/2) 
| vu G (we du = a + OP (46) 
0 


where I’ (x) is the Gamma function. The first term of (45) is obtained to be 





I(x, 2) = Cx-38-** sin ((k2 + ke) 2) (47) 
where 
dK" H, (0) A, k2 em /a {i+ ae 2k’2) (k? + k? ype 
war x H,(0) | VJ ok’ 


; 2 
+ MK k (2k, )Y8 tan (ke + ki) A) 


— As p,k’® (ki, ae k"2)1/2 (k? is kyr (2k° pr 
...(48) 
H, (0) = d, k* P (3/2) $i, (k’) . (49) 


H, (0) = (02k) (ke +R) (kO ~— 28) cos (KE + kt A)... (50) 


The scattered wave in (47) has the behaviour of a cylindrical wave close to the 
scatterer and a plane wave at distant points. 
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8. CONCLUSIONS 


The potential for the incident, reflected, Rayleigh and scattered waves in (3), 
(40)-(43) and (47) vanishes on the free surface (z = 0) of the ocean and the boundary 
condition (6) is satisfied. It can be further established that the condition (7) and all 
other boundary conditions are satisfied by the potentials of various waves. The hori- 
zontal component u of the displacement excited by various waves is absent on the free 
surface of the ocean. Thus the particle motion in the surface of the ocean is at right 
angles to it and the horizontal component of the displacement is absent. The waves 
reflected from the rigid plane boundary and then from the surface of the ocean are 
obtained in (41) having exactly half the amplitude of the incident wave. Multiply 
reflected generalised Rayeigh waves in (42) - (43) and the scattered waves in (47) are 
plane waves. When z is small, the scattered waves are of the form of exp (—k',r)[Vr,r 
being the distance from the corner of the plane boundary. Near the scatterered, the 
near-field behaves as a cylindrical wave. The scatted waves propagate with the speeds 
of the waves in the solid and not with the speed of the waves in the ocean. 


Numerical computations have been worked out for the scattered field and the 
generalised Rayleigh waves in the oceanic layer. For a basaltic ocean bottom, 


P. = 3P,, % = V3B,, B, = 3a: and h = 5.7km. 


Amplitude of the scattered wave is plotted versus the distance from the scatterer just 
below the free surface. It falls steadily as the waves proceed away from the scatterer 
(Fig. 3). The displacement excited by the scattered waves in the free surface falls rapidly 


65 
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55 
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20 





35 =40Km uF 


Fic. 3. Amplitude of the scattered wave just below the surface. 
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Fic, 4. Amplitude of the displacement excited Fic. 5. Amplitude of the displacement excited 
by scattered waves in the free surface. by Rayleigh waves versus depth. 
upto a distance of 10 km and then receds slowly when a = 1.52 km/sec, «2 = 7.95 


km/sec. and 8, = 4.56 km/sec. (Fig. 4). The amplitude of the displacement excited by 
Rayleigh waves is computed in terms of depth below the free surface and is found to 
fall down very slowly (Fig. 5). 
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A NOTE ON S-CLOSED SPACES 
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The class of s-closed spaces has been recently introduced by Di Maio and 
Noiri [Indian J. pure appl. Math., 18 (1987), 226-33]. In this paper it is 
pointed out that among weakly-7» spaces the class of s-closed spaces coincides 
with the class of S-closed spaces due to Thompson. We show that every 
infinite topological space is embeddable asa closed subspace in a connected 
S-closed space which is not s-closed. 


In a recent paper Di Maio and Noiri® have introduced and studied the class of 
s-closed spaces. They pointed out that s-closedness is equivalent to cd-compactness 
due to Carnahan? and to weak RS-compactness due to Hong’. The purpose of this 
note is to examine the relationship between the class of s-closed spaces and the more 
familiar class of S-closed spaces which was introduced by Thompson’ and studied by 
Cameron’. It is clear that every s-closed space is S-closed, but is the converse true ? 
No example of an S-closed space which is not s-closed has been provided in Di Maio 
and Noiri?, and according to Di Maio‘ this is an open question. In this paper we 
observe that these two classes of spaces coincide for a very large class of spaces, namely 
the weakly-7, spaces. We show that every infinite topological space can be represen- 
ted asa closed subspace of a connected S-closed space which is not s-closed. We 
also prove that there exist semi-7, S-closed spaces which fail to be s-closed. 


Let S be a subset of a space ¥. We denote the closure and the interior of S 
with respect to the space X by cly S and inty S respectively. A subset S of X is called 
semi-open (resp. regular closed] if S C clx (inty S) [resp. S = clx (intx S)]._ In parti- 
cular, every regular closed set is semi-open and the closure of every semi-open set is 
regular closed. The complement of a semi-open set is called a semi-closed set. The 
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semi-closure of a subset Sof X is the smallest semi-closed set containing S and is 
denoted by scly S. One easily verifies that scly S = S VU inty (clx S) for any subset 


S of X. 

Definition—A space X is called S-closed Thompson® [resp. s-closed®] if every 
semi-open cover of X contains a finite subfamily the closures [resp. the semi-closures] 
of whose members cover _X. 

It is obvious that every s-closed space is S-closed. Recall that a space X is said 
to be extremally disconnected if cly Uis open for each open subset U of X. It is 
observed in Maio and Noiri* that if X is exteremally disconnected then cly S = scly S 
for every semi-open set S in ¥. Thompson’ (Theorem 7) has proved that every Hausdorff 
S-closed space is extremally disconnected. Herrmann® (Theorem 3.7) has extended 
this result to the more general class of weakly-T» spaces Soundarajan®, i.e. spaces in 
which every singleton is the intersection of regular closed sets. Hence, among weakly-7, 
spaces the class of S-closed spaces coincides with the class of s-closed spaces. 


The question whether there is any difference between S-closed spaces and 
s-closed spaces is now settled by the following results. 


Theorem—Every infinite topological Space Y can be represented as a closed 
subspace of a space XY which is s-closed. 


Proor : Let Z be an infinite 7; space and let Z; = Z x {l} and Z, = Z x {2}. 
We may assume that Y ) (Z; U Z,) is empty. LetX¥ =Y UZ, U Z,. Fori = 1,2 
let Zi = {WiCZ,:W, = V x {i} for some open subset V of Z}. Let G73 = (Grex 
G=UUC, UC, where U C Y is open in Y and C; resp. C, are cofinite subsets of 
Z, resp. Z,}. It is easily checked that a topology on X is defined by taking BZ, U F, 
U Bs asa base. Clearly Z, and Z, are open in X and Y is a closed subspace of X. 
For i = 1, 2, we obviously have cly Z; = Y \) Z; and sclx Z; = Z;. For each y € be 
if Sy = Z, U {y} then Sy is semi-open in X and scly Sy = Sy. Hence {Sy: y € Y} 
U {Z,} is a semi-open cover of X which has no finite subfamily the semi-closures of 
whose members cover X¥. Thus Y is not s-closed. 


Corollary 1\—There exist S-closed spaces which are not s-closed. 


Proor : Let Z be an infinite S-closed T, space. Since Z; U Z, is dense in X, 
by Lemma 2.2 and Theorem 3.4 in Noiri’ it follows that X is S-closed. In particular, if 
Z is an infinite set Cainrying the cofinite topology then Y is even a connected S-closed 
T, space which fails to be s-closed. 


Recall that a space is said to be semi-T, if every pair of distinct points can be 
separated by disjoint semi-open sets. 


Corollary 2—There exist semi-T 2 S-closed spaces which are not s-closed. 


| PROOF : By our theorem and Corollary 1, if Y is an infinite discrete space and Z 
is an infinite Hausdorff S-closed space then XY is T; and S-closed but not s-closed. It 
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is easily checked that YX is semi-7>. For example, if »,¥, © Y and »:~y, then 


Z, U {y,} and Z2 U {y,} are disjoint semi-open sets containing y, and y, respectively. 
The other cases are treated in a similar fashion. 


ACKNOWLEDGEMENT 


We would like to thank the referee for valuable comments and suggestions. 


REFERENCES 


. Cameron, Proc. Am. Math. Soc. 72 (1978), 581-86. 
fe canbe, Ph. ~ Lae, University of Arkansas, 1973. 


x Herrmann, Proc, Am. Math. Soc. 75 (1979), 311-17. 

. C. Hong, J. Korean Math. Soc. 17 (1980), 39-43. 

. Noiri, Acta Math. Hungar. 35 (1980), 431-36. 

. Soundararajan, General Topology and its Relations to Modern Analysis and Algebra III (Proc. 
Penar Kampur 1968), Academia, Prague, 1971, pp. 301-306. 

9. T. Thompson, Proc. Am. Math. Soc. 60 (1976), 335-38. 





SUGGESTIONS TO CONTRIBUTORS 
The INDIAN JOURNAL OF Purg AND APPLIBD MATHEMATICS is devoted primarily to original 
research in pure and applied mathematics. 


< 
oa 


Manuscripts should be typewritten, double-spaced with sufficient - margins (including 
_abstracts, references, etc.) on one side of durable white paper. The initial page should contain the 
titel followed by-author’s name and full mailing address. The text should include only as much as 
is needed to provide a background for the particular material covered. Manuscripts should be 
submitted in triplicate. 


- 


é The author should provide a short abstract, in triplicate, not exceeding 250 words, summarizing 
the highlights of the principal findings covered in the paper and the scope of research. 


References should be cited in the text by the arabic numbers in superior. List of references 
should be arranged in the arabic numbers, author’s name, abbreviation of Journal, Volume number 
(Year) page number, as in the sample citation given below : 


For Periodicals * 
1. BR. H. Fox, Fund. Math. 34 (1947) 278. 


For Books’ . * : 
2. H.Rund, The Differential Geometry of Finsler Spaces, Springer-Verlag, Berlin, (1973) 
p. 283. ; 


Abbreviations for the titles of the periodicals should, in general, conform to the World List of 
Scientific Periodicals. 3 


All mathematical expressions should be written clearly including the distinction between capital 
and small letters. -Clear distinction between upper and lower cases of c,p,k,z,s, should be made 
while writing the expression in hand. Also distinguish between the letters such as “Oh’ and 
‘zero’; el) and 1 (one); v, V and v (Greek nu); r and y (Greek gamma); ¥, X and % (Greek chi); 
k, K and k (Greek kappa); Greek letter lambda (A) and symbol for vector product ( A); Greek 
letter epsilon (¢) and symbol for ‘is an element of? (€). The equation numbers are to be placed 
at the right-hand side of the page. The name of the Greek letter or symbol should be written in 
the margin the first time itis used. Superscripts and subscripts should be simple and should be 
placed accurately. ; 


Line drawings should be made with India ink on white drawing paper or tracing paper. _ 
Letterings should be clear and large. Photographic prints should be glossy with strong contrast. 
All illustrations must be numbered consecutively in the order in which they are mentioned in the 
text and should be referred to as Fiz. or Figs. Legends to figures should be typed ona separate 
sheet and attached at the end of the manuscript. 


Tables should be typed separately from the text and placed at the end of the manuscript. Table 
headings should be short but clearly descriptive. ; 


Proofs should be corrected immediately on receipt and returned to the Editor. Ifa large 
number of corrections are made in the proof, the author should pay towards composition charges. 
In case, the author desires to withdraw his paper, he should pay towards the composition charges, 


if the same is already done. 
For each paper, the authors will receive 50 reprints free of cost. Order for extra reprints 
- should be sent with corrected page proofs. 


Manuscripts, in triplicate, should be submitted to the Editor of Publications, Indian Journal of 
Pure and Applied Mathematics, Indian National Science Academy, Bahadur Shah Zafar Marg, 


New Delhi 110002 (India). 


RN 19797/70 | Regd. No. D-(C)-774. 


INDIAN JOURNAL OF PURE AND APPLIED MATHEMATICS 


No. 10 October 1988 ; Volume 19 
CONTENTS 
Page - 
Indefinite quadratic forms in many variables by MAry E. FLAHIVE Se 931 
On common fixed points in metric spaces by B. K. Ray Ss, re a OO 


On the stability of a system of differential equations with complex coeffi- 
cients by Z. ZAHREDDINE and E, F. ELsHEMAWEY aes ea 963 


Singularly perturbed initial value problems for differential equations ina 


Banach space by N. RAMANUJAM and V. M. SONANDAKUMARI ad 973 


Sums involving'the largest prime divisor of an integer II by JEAN-MAriE Deg 
KONINCK and R. SITARAMACHANDRARAO we ay Ke 990 


On the problem of three gravitating triaxial rigid bodies by S. M. / 
ELSHABOURY iA vie re sn (3 a 1005 


On a particular initial value problem, with an application in reservoir 


analysis by B. HOFMANN an $ sei aaa san = erp OeE 
P-Wave scattering at a coastal region in a shallow ocean by P. S. DesHwaL ; 

and NARINDER MOHAN Sad re yy .? orate 1020 
A note on S closed spaces by MAXIMILIAN GANSTER and Ivan L. REILLY ... 1031 


Published and Printed by Dr O. N. Kaul, Executive Secretary, Indian National Science Academy, 
Bahadur Shah Zafar Marg, New Delhi 110002, at Leipzig Press, D-52,;N.D.S.EB. Part I 
New Delhi 110049, Ph. 622490 ; 


Caste 


